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FOREWORD 


This  report  describes  work  accomplished  in  the  program  "Quasi 
Area  Rule  for  Heat  Addition  in  Transonic  and  Supersonic  Flight  Regimes" 
conducted  under  USAF  MIPR  APO-71-007.   The  work  was  accomplished 
during  the  period  1  Oct  1970  through  31  Aug  1971  at  the  Naval  Post- 
graduate School,  Monterey,  California,  under  the  direction  of  Dr.  Allen 
E.  Fuhs.   The  report  was  submitted  on  1  Nov  1971. 

The  program  was  sponsored  by  the  Air  Force  Aero  Propulsion  Labora- 
tory, Wright-Patterson  Air  Force  Base,  Ohio  under  Project  3066,  Tur- 
bine Engine  Propulsion,  Task  306603,  Advanced  Component  Research. 
Dr.  Kervyn  D.  Mach,  AFAPL/TBC,  Turbine  Engine  Components  Branch, 
was  the  project  engineer. 

Publication  of  this  report  does  not  constitute  Air  Force  ap- 
proval of  the  report's  findings  or  conclusions.   It  is  published 
only  for  the  exchange  and  stimulation  of  ideas. 


ERNEST  C.  SIMPSON     *~~ 
Director,  Turbine  Engine  Division 
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ABSTRACT 

Body  shapes,  including  axisymmetric  and  three  dimensional,  have  been 
developed  to  minimize  wave  drag.   The  von  Karman  ogive  and  the  area  rule  are 
examples.   Similar  work  has  not  been  accomplished  for  optimum  shapes  with 
propulsion.   Propulsion  can  be  divided  into  two  categories — those  devices  with 
internal  heat  addition  and  those  with  external  burning.  For  internal  heat 
addition  an  analytical  model  is  formulated  which  introduces  the  propulsive  disc, 
Attention  is  shifted  to  external  burning,  which  is  examined  for  one  dimensional 
and  two  dimensional  linearized  flow.   Heat  fronts  and  combustion  fans  are 
discussed  as  examples.  Forces  on  a  heat  source  in  a  uniform  stream  and 
adjacent  to  bodies  are  derived.  Several  possible  applications  are  examined 
including  base  pressure  augmentation  by  external  burning,  spin  recovery  using 
external  burning,  and  transonic  boattail  drag  alleviation.   Previous  work  on 
base  pressure  augmentation  has  used  a  two  dimensional  planar  model.  A  two 
dimensional  axisymmetric  model  is  examined. 
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QUASI  AREA  RULE 

FOR  HEAT  ADDITION  IN  TRANSONIC 

AND  SUPERSONIC  FLIGHT  REGIMES 

Dr.  Allen  E.  Funs 
Professor  of  Aeronautics 
Naval  Postgraduate  School 

INTRODUCTION 

In  subsonic  flow,  drag  can  be  attributed  to  viscous  forces  on  the  body- 
surface  and  to  the  release  of  vorticity  in  the  wake.   The  vorticity  is  associated 
with  the  generation  of  lift.   In  supersonic  flow,  there  is  a  new  phenomenon  lead- 
ing to  drag;  it  is  the  radiation  of  waves  from  the  body.   These  waves  carry 
energy  whose  origin  is  the  work  done  by  drag  and  vehicle  motion. 

For  almost  every  case  of  a  body  flying  through  a  medium,  drag  is  considered 
to  be  a  degradation  of  performance.  Consequently,  considerable  effort  has  been 
devoted  to  optimum  shapes  giving  minimum  drag. 

One  approach  for  optimizing  a  body  shape  is  to  represent  the  body  by  sources 
and  sinks.     For  a  closed  axisymmetric  body,  all  sources  and  sinks  are  on  the 
axis.   It  is  found  that  for  a  nonlifting  body,  the  slope  of  the  cross  sectional 
area  curve  equals  the  distribution  of  sources  and  sinks,  f (x).   Distance  along 
the  axis  is  x.   By  means  of  the  calculus  of  variations,  one  can  arrive  at  shapes 
giving  minimum  wave  drag.   There  are  constraints  on  the  body  geometry,  and  the 
optimum  shape  depends  on  these  constraints.   Two  examples  are  given  in  Table  I. 
The  shapes  given  in  Table  I  are  for  minimum  wave  drag.  Viscous  effects  are  not 
included,  and  the  bodies  are  nonlifting. 


Table  I.   Optimum  Shapes 


Name  of  Body  Constraints  Wave  Drag 


Von  Karman*  Ogive ^  '     Given  length,  l>    and  base  area,  S       =qf  —  ] 

At2/ 

Sears**Ogive^'         Given  length,  lf    and  volume,  V         !££q  [  X_ 

*Used  on  Low  Altitude  Short  Range  Missile. 
**Discovered  by  W.  R.  Sears  of  Cornell  University. 


It  was  discovered  by  Lomax  and  Heaslet    and  by  Whitcomb    that  the 
optimum  shape  for  an  axisymmetric  body  applies  to  slender  bodies  with  wings  and 
rudders  if  one  interprets  the  area  distribution  properly.   It  is  the  area  cut  by 
the  Mach  cone.   If  one  uses  this  area  distribution  from  a  three-dimensional  body 
to  match  the  area  distribution  of  an  optimum  axisymmetric  body,  then  minimum  wave 
drag  occurs.   This  procedure  is  known  as  the  area  rule. 

For  a  body  with  thrust,  there  are  several  new  variables  introduced.  Thrust 
may  be  developed  by  adding  heat  to  a  flow  internal  to  the  body.  Thrust  also  may 
be  obtained  by  adding  heat  external  to  the  body. 

Consider  external  burning.   Is  there  an  area  rule  for  a  body  with  external 
heat  addition  analogous  to  the  area  rule  involving  body  shape  only?  For  this  case, 
there  are  two  functions  to  optimize.   One  is  body  shape,  and  the  other  is  the 
distribution  of  heat  addition.   Table  I  gives  some  typical  constraints;  an 
additional  constraint  for  the  optimization  could  be  specific  fuel  consumption. 


Instead  of  finding  minimum  drag,  the  procedure  should  give  the  best  thrust- 
minus-drag  or  zero  thrust-minus -drag  depending  on  the  constraints.   Since  both  the 
distribution  of  heat  and  cross  sectional  area  are  involved,  an  appropriate  name 
is  "quasi  area  rule." 

When  thrust  is  developed  by  adding  heat  to  a  flow  internal  to  the  body,  the 
quasi  area  rule  suggested  in  the  previous  paragraphs  does  not  apply.   For  in- 
ternal burning,  the  waves  due  to  heat  addition  obviously  do  not  occur  as  in  the 
case  of  external  burning.   There  is,  however,  a  change  in  stagnation  pressure 
and  temperature  between  the  inlet  and  nozzle.   Also  the  velocity  vector  changes 
from  inlet  to  exhaust.   An  energy  disc  can  represent  those  changes.   This 
report  discusses  aspects  of  the  thrust-minus -drag  problem  for  an  airbreathing 
engine  with  confined  heat  addition.   Figure  1  illustrates  the  various  problems 
that  can  be  studied  with  internal  burning  when  the  emphasis  is  on  thrust-minus- 
drag.   The  analysis  is  linearized.   Two-dimensional  (2D)  planar  is  the  geometry 
of  a  propulsive  wing.   Axisymmetric  2D  geometry  is  an  isolated  nacelle  without 
fins  or  wings.   Addition  of  fins  or  wings  yields  the  2D  propulsion  with  3D  body. 
The  connection  between  these  two  geometries  is  the  area  rule. 

Figure  2  outlines  the  problems  of  interest  for  external  burning.   The 
analysis  is  linearized.  For  the  2D  planar  case,  the  heat  and  mass  sources  can 
be  distributed  over  a  surface  or  throughout  a  volume.  Volume  distributions  must 
meet  requirements  for  slender  body  theory.   Axisymmetric  2D  bodies  are  related 
to  3D  bodies  by  the  quasi  area  rule. 

The  linearized  2D  planar  flow  with  heat  addition  has  been  studied  extensively. 
A  solution  to  the  2D  axisymmetric  case  does  not  appear  in  the  scientific 
literature.   Solutions  for  3D  bodies  with  heat  addition  have  not,  as  yet,  been 
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obtained.   In  order  to  develop  the  quasi  area  rule,  these  solutions  are 
essential  building  blocks. 

PROPULSION  WITH  INTERNAL  HEAT  ADDITION 
Introduction 

For  a  supersonic  aircraft,  the  variation  in  nozzle  pressure  ratio  (NPR) 
is  very  large.   This  requires  a  variable  geometry  nozzle.   Under  certain  flight 
conditions,  the  nozzle  setting  and  nozzle  flow  result  in  large  base  or  boattail 
drag. 

One  motivation  for  treating  the  propulsion-airframe  problem  as  a  unit  from 
the  area  rule  point  of  view  is  to  obtain  new  design  approaches.   First,  the 
typical  variation  in  NPR  and  nozzle  geometry  will  be  calculated  for  a 
Mach  2.5  aircraft.   These  results  will  be  discussed  to  illustrate  the  problem. 
Second,  the  model  for  representation  of  a  propulsion  system  with  internal  heat 
addition  will  be  developed.   Third,  approaches  to  solution  of  the  model  will  be 
outlined. 
Exhaust  Variation  with  Mach  Number 

The  aim  of  this  section  is  to  calculate  nozzle  pressure  ratio,  nozzle  area 
ratio,  mass  flow  ratio,  and  thrust  of  an  ideal  turbojet  as  a  function  of  Mach 
number.   Knowledge  of  the  trends  is  essential  in  understanding  aircraft  exhaust 
problems. 

Consider  an  ideal  turbojet  with  an  ideal  inlet.   There  is,  however,  a 
restriction  on  the  Mach  number,  M?,  into  the  compressor.   The  thrust  is  given  by 

F  =  f—  1^  )fc,g  (1) 

\V]\md*l    d 

where  m  is  a  mass  flow  of  air,  and  fh,  is  mass  flow  of  inlet  at  design  conditions, 
a  'a 


For  the  inlet,  the  various  areas  are  illustrated  in  Figure  3(a).   Station  0 
is  at  a  point  in  the  freestream  where  pn  is  ambient  pressure.   Stations  1  and  2 
are  inlet  throat  and  compressor  face  respectively.   The  mass  flow  ratio  is 

Ad "  Moa  At  H 

where  Mn  is  flight  Mach  number,  and  M_,  is  design  flight  Mach  number.   Areas 
An,  A  ,  and  A§   are  identified  in  Figure  3(a).   Both  a  fixed  area  and  variable 
geometry  inlet  are  considered.   The  ratio  A.VA  is  obtained  from  isentropic 
flow  tables  for  given  M_  for  the  variable  geometry  inlet.  For  a  fixed  geometry 
inlet  A-j/A.  comes  from  flow  tables,  entering  the  flow  tables  with  the  Mach 
number  downstream  of  a  normal  shock  at  M-.   The  mass  flow  ratios  are  shown  in 
Figure  4  for  both  cases.   A  variable  geometry  operates  supercritically  and 
swallows  more  air.  A  variable  geometry  ramp  or  spike  inlet  may  have  superior 
mass  handling  capability  at  subsonic  speeds. 

The  thrust  per  pound  of  air  for  an  ideal  turbojet  is  given  by 
F 


wa 


rv-^  -*r--i} 


-1         (3) 


where  the  symbols  not  previously  defined  have  the  following  meaning: 
an  =  freestream  speed  of  sound 

T  =  TTn/Tn  =  ratio  of  stagnation  to  static  freestream  temperature 
T  =   TT,/TTp  =  stagnation  temperature  ratio  across  the  compressor 

The  symbol  7**  deserves  more  discussion.   It  is 

r*-   rrrcrb=  VT0  M 

where  T,    =  T^./T--  =  stagnation  temperature  ratio  across  the  burner  or  combustor. 

The  value  of  TT/  is  limited  and  has  a  maximum  value.   Hence  for  a  given  engine 

and  altitude,  T*   is  fixed  when  the  engine  is  at  maximum  power  or  maximum  TIT. 
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FIGURE  3.  NOTATION  FOR  TURBOJET  CALCULATIONS. 
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FIGURE  4.  MASS  FLOW  OF  AIR  AS  A  FUNCTION  OF  FLIGHT  MACH  NUMBER 


The  station  numbering  is  shown  in  Figure  3(b).   The  throat  area  for  the 
nozzle  is 

I   PT8 

neglecting  fuel  addition.   R  is  the  gas  constant.  T  is  a  function  of  the  ratio  of 
heat  capacities,  Y»  and  has  a  value  of  .683  for  y  =  1«4.«   The  nozzle  pressure  ratio, 
NPR,  is  pTo/pn.   A  ratio  of  A~  to  the  design  value  is 


A-     h       NFR   /T 

-2L=  J  — ^i/—  (6) 

A?d    ^  NFH  VTTgd 

The  area  Aft  is  obtained  from  NPR,  isentropic  flow  tables,  and  the  assumption 

Pg  =  Pq.   Two  useful  ratios  are  A«/A  and  A_/A  ,  i.e.,  comparing  nozzle  areas  with 

compressor  face  area. 

Calculations  were  made  for  the  values  shown  in  Table  II. 

One  of  the  results  of  the  calculation  is  NPR  which  is  shown  in  Figure  5.   Some 
of  the  other  results  are  summarized  in  Table  III   The  symbols  for  the  turbine  ratios 

are  Tt  =  TT5/TT4  and  nt  =  pT5/pTr 

Figure  6  illustrates  the  area  variations  required  for  a  supersonic  turbojet 

operating  at  maximum  TIT  for  different  K^.   Examination  of  Figure  6  and  the  data  in 

Table  III  shows  why  turbojet  nozzles  usually  operate  underexpanded  at  design  Mach 

number.   The  percentage  variation  in  A,  is  38  per  cent  based  on  design  Ay.   The 

percentage  variation  in  A«  is  77  per  cent  based  on  design  Ag.   NPR  varies  from  U   at 

low  subsonic  flight  to  more  than  25  at  design  M~  0f  2.5. 
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Table  II.   Values  for  Calculation  of  Turbojet  Quantities 


Maximum  turbine  inlet  temperature 
Altitude 

Frees tream  density 
Freestream  speed  of  sound 
Compressor  area 
Tip  radius  at  face 
Hub  radius  at  face 
Maximum  Mach  number  at  compressor 
Inlet  throat  area  at  design 
Design  flight  Mach  number 
Compressor  pressure  ratio  at 
100  per  cent  corrected  speed 


U 


'0 


M 


2max 


ltd 


M. 


'0d 


2000°F 


Sea  Level 


.002378  slugs/ft 

1117  ft/sec 

8.0$  ft2 

2  ft 

1.2  ft 

0.65 

7.1  ft2 

2.5 


8 


11 


25 


20 
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FIGURE  5.  NOZZLE  PRESSURE  RATIO  AS  A  FUNCTION  OF  MACH  NUMBER. 
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Table   III. 

Results 

of  Calculat; 

Lon  for 

Turbojet  Perfo: 

nuance  and 

Areas 

Mo 

NFR, 
a 

NPR 

A7 
*7d 

h 

A 
c 

A8 

A 
c 

Tt 

\ 

0.2 

5.98 

.625 

.421 

.526 

.829 

.518 

0.4 

5.63 

.627 

.423 

.542 

.825 

.508 

0.6 

5.05 

.628 

.425 

.574 

.818 

.495 

0.8 

4.42 

.637 

.427 

.620 

.808 

.475 

1.0 

3.77 

.650 

.434 

.677 

.796 

.45 

1.2 

3.15 

.663 

.443 

.762 

.780 

.42 

1.4 

2.58 

.702 

.448 

.860 

.765 

.391 

1.6 

2.14 

.712 

.477 

1.02 

.742 

.352 

1.8 

1.77 

.750 

.503 

1.20 

.719 

.315 

2.0 

1.47 

.792 

.537 

1.42 

.694 

.278 

2.2 

1.24 

.861 

.582 

1.72 

.665 

.240 

2.4 

1.07 

.950 

.642 

2.07 

.634 

.203 

2.5 

1.00 

1.00 

.685 

2.28 

.616 

.183 
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As  shown  in  Figure  6,  to  obtain  ideal  expansion  the  nozzle  is  converging- 
diverging  for  all  Mn  illustrated.   Usually  the  nozzle  is  a  simple  converging 
nozzle  from  subsonic  through  transonic  range  of  ML.   This  results  in  a  large 
boattail  area.   At  the  design  M„  the  exit  area  for  ideal  expansion  significantly 
exceeds  the  maximum  area  of  the  propulsion  system.   The  nozzle  is  cut  off,  relative 
to  ideal,  at  some  area  appropriate  to  maximum  engine  area.   The  result  is  under- 
expansion  and  a  ballooning  plume.   A  shock  wave  may  originate  at  the  diverging 
streamlines  of  the  plume,  giving  a  drag  increment. 
Model  for  Propulsion  with  Internal  Heat  Addition 

Consider  a  propulsion  device  as  sketched  in  the  top  of  Figure  7.   To 
represent  the  wave  drag  of  the  inlet,  it  is  necessary  to  duplicate  streamline 
A-A.   It  does  not  matter  what  the  details  are  of  the  inlet  geometry  internal 
to  the  streamtube  bounded  by  A-A. 

The  inner  body  which  generates  a  flow  having  streamline  A-A  is  shown.   Note 
that  it  is  not  the  same  shape  as  the  inlet  spike.   A  distribution  of  sources  on 
the  axis  could  generate  the  inner  body  as  well  as  streamline  A-A.   Note  that  the 
spike  shock  wave  does  not  have  the  same  angle  as  the  cowl  shock  wave. 

If  one  applies  the  momentum  relations  to  the  control  volume  of  Figure  7, 
the  thrust  minus  drag  is  given  by 


T  -  D       1       I  1 

A 


(u2  +  u2)dA 
y         r' 


DR 

(p„  -  p    )dA  -  -2  +  2      luudA 


(7) 
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FIGURE  7  MODEL  OF  AN  INLET  USING  SOURCES 
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where  u  $   u  ,  and  u  are  the  dimensionless  perturbation  velocities  arising 
x   r       y  v  6 

from  the  nacelle.   It  is  apparent  that  for  the  case  where  the  exhaust  exit 

defines  a  plane,  as  illustrated  in  Figure  7,  the  thrust  and  drag  can  be 

clearly  separated.   Waves  generated  by  the  exhaust  plume  do  not  influence  T  -  D. 

Following  the  analysis  suggested  by  the  model  of  Figure  7,  one  can  obtain 
a  nacelle  of  known  wave  drag.   Figure  8  illustrates  the  concept.   Start  with  an 
ogive  of  known  shape.   Flot  the  streamlines  in  the  flow  adjacent  to  the  ogive. 
Fick  the  radii  r~  and  rR  suitable  for  the  nacelle.   Insert  the  nacelle  within 
the  streamtube.   The  wave  drag  of  the  ogive  and  nacelle  will  be  the  same.   The 
drag  can  be  calculated  using  the  potential  function  for  the  ogive.   In  the 
lower  part  of  the  figure,  a  nacelle  with  a  spike  inlet  is  shown. 

Consider  some  thrust  device  inside  a  momentum  control  volume  shown  in 
Figure  9.   How  does  one  modify  the  variables  so  as  to  change  conditions  on  the 
surface  of  control  volume  and  improve  T  -  D?  What  does  one  have  available  to 
change  the  wave  drag?  Available  variables  are:   (a)  number  and  location  of 
nozzles,  (b)  velocity  profile  at  nozzle  exit,  and  (c)  body  geometry. 

To  obtain  an  optimum  configuration,  one  must  have  a  method  to  represent 
the  body  and  exhaust  analytically.   The  analytical  representation  can  be 
subjected  to  the  calculus  of  variations. 

The  body  can  be  represented  by  sources  and  sinks.   Changes  in  stagnation 
pressure,  stagnation  temperature,  velocity  vector,  and  other  flow  quantities 
make  representation  of  exhaust  more  difficult. 

Carrying  the  idea  of  Figure  7  further,  let  us  look  at  an  ideal  ramjet. 
For  an  ideal  ramjet,  the  following  equalities  hold  true: 
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p6 

pr 

T 
Ti 

TT6 

pi  = 

p6 

*1"  Ao 

PT6  =  Pn 

VC  "  TT1 
For  an  ideal  ramjet  (or  for  any  other  propulsion  device  satisfying  the  above  set 
of  equalities),  one  can  represent  the  flow  by  sources,  sinks,  and  an  "energy 
disc."  This  is  illustrated  in  Figure  10;  stations  1  and  6  are  identified  in 
Figure  10.   The  inner  body  duplicates  streamline  ABC.   The  inner  body  is  formed 
from  the  source-sink  distribution  in  the  right-hand  column  of  Figure  10.   The 
change  in  density,  velocity,  and  stagnation  temperature  is  achieved  by  an 
"energy  disc."   In  fact,  that  defines  the  energy  disc.   It  changes  the  flow  so  that 
the  equalities  list  above  is  correct.   The  sources,  sinks,  and  energy  disc  shown 
in  the  bottom  of  the  right  column  would  have  the  same  thrust  and  drag  as  the 
ideal  ramjet  sketched  in  the  left-hand  column. 

This  model  looks  as  if  it  is  satisfactory  for  an  ideal  ramjet,  but  it  has 
some  drawbacks: 

1.  The  energy  disc  is  not  a  heat  front  or  other  simple  fluid  flow. 

2.  When  A.  /  A/,  it  may  be  difficult  to  apply  the  energy  disc. 

3.  For  the  case  p_,,  j-   p_,  ,  some  means  of  changing  p_,  must  be  added. 

The  energy  disc  is  bounded  by  a  streamtube.   Changing  pT  within  the  disc  gives  a 
mismatch  across  the  disc  in  regard  to  streamline  shape. 

A  satisfactory  way  of  specifying  the  conditions  at  station  6  is  needed,  and 
the  energy  disc  has  been  suggested.   One  wants  to  retain  sources  and  sinks  for 
two  reasons.   Wave  drag  is  readily  calculated  if  you  know  the  distribution  function. 
The  distribution  function  can  be  optimized  to  yield  minimum  wave  drag. 
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The  energy  disc  represents  the  propulsion  device.   A  particular  engine 
generates  changes  in  p-,,  T„,,  and  velocity  between  the  inlet  and  exhaust.   For 
the  case  of  a  plane  for  an  exit  as  shown  at  station  6  in  Figure  10y  one  does  not 
need  the  complexity  of  an  energy  disc.   Equation  (7)  indicates  that  drag  and 
thrust  are  separable  for  this  case.   When  one  varies  quantities  (a),  (b),  and 
(c)  discussed  previously  as  available,  independent  variables,  then  the  energy 
disc  may  be  an  attractive  model. 

This  section  is  closed  by  noting  that  a  model  has  been  proposed  but 
has  not  been  developed. 

EXTERNAL  BURNING;  BASIC  EQUATIONS  AND  SOLUTIONS 
One-Dimensional  Heat  Addition 

To  provide  a  basis  for  comparison  with  subsequent  sections,  one-dimensional 
heat  addition  will  be  briefly  reviewed. 

There  are  several  approaches  to  one-dimensional  heat  addition.   For  unsteady 

(6) 
flow  there  is  the  self  similar  technique  developed  by  Sedov.     We  consider 

steady  flow  only  here.   There  is  the  famous  solution  for  constant  area  heat 

addition,  i.e.,  the  solution  by  Rayleigh.   Rayleigh  flow  is  discussed  in  most 

(7) 

gas  dynamics  books.     When  there  is  simultaneous  variation  of  area  and  heat 

addition,  the  techniques  outlined  in  the  Princeton  Series,  Volume  III,  are 

(8) 
applicable. 

(9) 
The  influence  coefficients  developed  by  Shapiro  and  Hawthorne    will  be 

discussed  here.   Table  IV  summarizes  the  influence  coefficients  for  heat 

addition  and  isentropic  flow  with  area  change.   The  influence  functions  relate  a 

fractional  change  of  a  dependent  variable,  say  M  ,  to  a  fractional  change  of  the 

independent  variable,  A  or  h_,.   So  far  as  heat  addition  is  concerned,  the 
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Table  IV.   Influence  Coefficients  for  Area  Change  and  Heat  Addition 


dhT 

p 

dM2  -2(1  +  ^-=-^M2)  1  +  yM2 

du 

u 

dT 
T 

d£ 
P 


dj> 

p  1  -  M2  1  -  M2 


dF 
F 

ds 
c 


dA 

T 

-2(1    +  ^lM2) 

2 

1-M 

1 

1   -  M2 

(y  -  Dm2 

1-M2 

M2 

1-M2 

yM 

1-M2 

1 

1  ♦  ym2 

1-M2 

1 

2 
1  -  IT 

l  -  yM2 

1-M2 

1 

'l-M2 

yM2 

0 


velocity,  u,  and  static  pressure,  p,  change  as  they  would  flow  into  decreas- 
ing area.  For  u  and  p,  heat  addition  and  area  variation  produce  the  same 
result. 

For  unconfined  flow  with  heat  addition,  there  are  some  analogies  between 
the  flow  due  to  a  particular  body  shape  and  the  flow  due  to  a  given  distribution 
of  heat  addition.   These  analogies,  which  are  limited  in  scope,  will  be  pointed 
out  as  they  arise  in  subsequent  discussion.   The  fact  that  heat  addition  and  a 
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body  may  influence  the  flow  the  same  way  may  be  understood  using  the  du/u  and  dp/p 
influence  coefficients. 

The  influence  coefficients  are  the  first  terms  in  a  Taylor  series  expansion 
of  the  dependent  variables  in  terms  of  the  independent  variables.   Many  of  the 
results  of  linearized  heat  addition  in  an  unconfined  flow  can  be  related  to 
influence  coefficients.   This  statement  is  applicable,  for  examble,  to  heat 
fronts  and  combustion  fans. 
Linearized  Equations  for  Combined  Mass  and  Heat  Addition 

The  equations  of  motion  for  steady  flow  with  heat  and  mass  addition  are  as 
follows: 

Continuity  div(pu )  =  m  (8) 

Momentum  p—  +  vp+mu=0  (9) 

_^     2 

Energy  div  pif(h  +  j-)  =  Q  (10) 

In  the  preceding  equations,  m  is  the  rate  at  which  mass  is  added  per  unit 
volume.   Likewise  Q  is  the  rate  at  which  heat  is  added  per  unit  volume. 

To  linearize  these  equations,  one  introduces  perturbation  quantities  which 
are  defined  as  follows: 

u=e*(U  +u,)+"e"v'  +ew' 
x  y     z 

When  the  perturbation  quantities  are  introduced,  the  resulting  equations  are 

Continuity  p  u1  ♦  p^v'  +  Up '  =  m  (12) 

co   x    °°   y     x 

P«,u»uf  +  p1  =  -  mU_ 
/  °°   x    x  t      \ 

Momentum  -j  (13; 

P^UMv'  ♦  p»  =  0 
roo  «o  x    y 


2A 


P„U«  8  /P'   YP'\    am   2 

Energy      -=—=■  t- ♦  -^-( r  -  VT     =  Q  H 

We  have  formulated  a  2D  planar  flow  so  that  w'  =  0.   To  avoid  writing 
subscript  oo,  we  now  drop  it  in  subsequent  equations.   The  notation  with 
subscript  x  means  partial  differentiation  with  respect  to  x» 

Starting  with  Equations  (8)  to  (10),  one  could  obtain  the  linearized 
equations  appropriate  for  heat  and  mass  addition  in  axisymmetric  flow.   The 
equations  have  been  solved  for  mass  addition  (sources)  but  not  for  heat 
addition  in  axisymmetric  flow. 

The  solution  for  a  source  in  two-dimensional  planar  flow  will  now  be 
obtained.   Consider  a  row  of  point  sources  evenly  distributed  along  the 
z-axis.   Axis  orientation  relative  to  the  flow  is  shown  in  Figure  11.   For 
a  point  F(x  -  x,  ,  y  -  y,  )>  only  the  sources  within  the  Mach  forecone  influence 
point  P  as  shown  in  Figure  11.   To  find  the  potential  in  two-dimensional  flow, 
add  up  the  sources  along  the  z-axis  from  k  to  d.   Points  k  and  d  are 
illustrated  in  Figure  11.   From  the  geometry  of  Figure  11 

zd  =  (x  -  xx)  tan  u  -  (y  -  y^  (15) 

Integration  of  a  row  of  point  sources  yields   the  potential 

0(x,y,O)  =   -  (16) 


-  z 


2Try(x  -x1)2  -  p2^-  p2(y  -ya) 
a 

This   integrates  to 
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FIGURE  II.  POTENTIAL  FOR  A  LINE  SOURCE  IN  SUPERSONIC  FLOW. 
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1     _x  pZj 

0(x, z)  =  -  —  sin 


zd 


(17) 


P"        -y(x  -   xx)  -  £2(y  -  y1)2 

Inserting  the  limits  gives 

0(x,y)  =  ^  (18) 

The  potential  is  a  constant.   Think  of  a  series  of  Mach  cones  with  apex 

at  the  z-axis.   Upstream  of  these  cones,  the  potential  is  zero.   Downstream 

of  the  envelope  of  Mach  cones,  it  is  constant.   One  can  represent  this  by 

0(x,y)  =  -  |j  I  [(x  -  xx)  1  p(y  -  yx)]  (19) 

where  I  is  the  unit  function.   The  velocity  components  are,  for  a  unit 
source  strength, 

u  ■  H  ■  -  ksh  -  *i>  ±  p(y "  yi}]  (20) 

and 

w  =  Sf =  7 26[(x  -  xi)  -  ?(y  -  ^  <21> 

The    (  +  )   sign   in  front  of   f  gives   right-running  waves,   whereas   the    (-)   sign 
gives   left-running  waves. 

The   solutions  for  a  heat  source   in  supersonic  flow  are 

U'  =    -   (2Yp/}q  6(x  "  *y)  (22) 

V  =    {r^plhS(x  -  fr)  (23) 

p,,   (T-al?Mq    6(x.  py)  (at) 

pl.    (Y  -  l)Mq    6(x  .  py)   -  (T  :1)q6(y)l(x)  (25) 

2a  p  a  U 
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One  can  verify  that  there  are  the  appropriate  solutions  when  m  =  0  and  Q  /  0 
by  substitution  into  Equations  (12)  to  (14-).   The  second  term  in  Equation  (25) 
is  a  wake.   Pressure  in  the  wake  is  the  ambient  value,  and  the  velocity 
returns  to  freestream  value.   The  solutions  in  Equations  (22)  to  (25)  are 
valid  for  all  x  and  y  including  the  source  point  x  =  y  =  0.   The  symbol  q, 
which  has  units  heat/(length)(time),  will  be  discussed  shortly.   The  line  heat 
source  in  a  planar  flow  is  illustrated  in  Figure  12. 

Similar  solutions  can  be  written  for  the  case  of  mass  addition;  these  are: 

u«  =  -  i  6(x  -  py)  (26) 

v«  =  ~6(x-py)  (27) 

P'  =  §f&(x-Py)  (28) 

P'  =  ^6(x-^y)  (29) 

Substitution  of  Equations  (26)  to  (29)  into  Equations  (12)  to  (H) 

indicate  these  are  solutions  except  at  the  mass  source  point  x  =  y  =  0. 

With  Equations  (22)  through  (29),  one  has  solution  for  volume  distribution 

of  heat  and  in  regions  where  m  =  0.   This  is  adequate  to  describe  planar 

bodies  with  adjacent  heat  addition  zones. 

Volume  Heat  Addition 

Since  the  heat  may  be  added  throughout  a  volume  rather  than  a  line  source 

(the  z-axis  in  Figure  11),  it  is  worthwhile  to  formulate  heat  as  a  volume 

release.   Define,  as  in  Equations  (10)  and  (14.) 

D  _  heat 

(span) (area) (sec) 
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3$ 


line   ^^ 
heat  source 


STREAMLINE 


wake  with  density 
P=P<o-P, 


Mach  angle /x=arcsin  M 


infinitesimal  width 


FIGURE  12.  SMALL  PERTURBATION  HEAT  SOURCE  IN  A  SUPERSONIC 
FLOW. 
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Then  we  define 

heat 
^   (span) (sec) 

q  =  Qdxdy 

Another  area  element  may  be  better  than  dxdy.   Some  element  with  its  sides 

parallel  to  a  Mach  line  is  preferred.   We  think  of  Q  =  Q(x,y).   Figure  13 

is  a  series  of  drawings  illustrating  the  variables. 

To  get  heat  per  unit  span 

L  rl 


H  =    I    I   Qdxdy  =  QL*  (^0) 

*   0  J   0      L if  Q  =  constant 

Heat  released  in  the  complete  volume  is  U  =  Hs. 

~.  _  heat       „  _     heat 

sec  (span) (time) 

Now  to  tie  into  the  2D  flow  equations,  we  use  Equation  (22)  as  an  example, 

»'  ■  -  (Ta^)q  &(x  -  rr)  fa) 

At  a  point  on  a  characteristic,  u*  results  from  a  source  on  that  same 
characteristic.   If  there  is  more  than  one  source  on  that  characteristic, 
one  adds  the  influence  of  each. 

W  -   lXq,  ♦  q2) 
«'■  ^ ^6(x-py)  (12) 

For  many  heat  sources,  in  fact  n  heat  sources, 

i  =  n^ 

i  =  1 
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^heat  soucef^ 
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■>  U 


(E) 


point  P 


FIGURE   13. 
L LUSTRATION  OF  QUANTITIES  RELATED  TO  SPECIFICATION  OF  HEAT  RELEASE. 


The  equation  for  u1  is  correct  if  one  thinks  of 

_  Q,    heat 

^i     ^   (span ) (length) (sec ) 

when  the  heat  source  is  distributed  in  the  volume.   When  heat  sources  are  not 

discrete  but  are  distributed,  then  one  wants  to  integrate  along  a  characteristic. 

Let  s  be  distance  along  the  characteristic.  We  express  q.  in  terms  of  Q  and  ds 

q.  =  Q  sin  u  ds  =  dq  (34.) 

The  perturbation  velocity  is 


r3 

%$   1     «^M--^ 


Q(x,y)ds  (35) 


u'  =  - 

0  "  0 

Since  s,x,y  are  related  quite  simply 

x  =  s  cos  ji,      y  =  s  sin  p.,  x  =  fy 

one  can  integrate  u1  several  ways: 

fs  C 

sin  u  J   Q  s  cos  ju,  s  sin  p.   ds    or   j   Q(fy,y)dy 
0  0 

f 

or      Q(x,x/p)dx/e 

Referring  to  the  region  bounded  by  0  <  x  <  L  and  0  <  y  <  I   in  Figure  13(b) 
and  for  case  Q  =  constant  inside  the  region 


(36) 
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X 

f 


sin  u  Qds  =  Ql  =   I  Qdy  =  Ql  = 
0  0  "0 

The  pressure  perturbation  would  be 


vdx 


Q(x,x/f)y  =    (tanja)pQl  =    <tf  (37) 


.  .  (r   -l)MQt 
p       2a  p 

In  terms  of  Q(x,y),  the  perturbation  quantities  are 


(38) 


...  ,(Y-D 

2rPg 


Q(x,y)  sin  ju  ds 


(39) 


s 
1  =  ^p"^  I  Q(x»y)  sinju  ds 


Uo) 


s 


(u) 


p'  = 


2a3p 


Y  -  1 

Q(x,y)  sin  u  ds  -  2— 5 — 

a  U 


Q(x«,y')6(y  -  y')X 


x'y' 


I(x  -  x»)dx'dy' 
The  wake  is  given  by  the  integral  above;  as  an  example  for  x  >  x.  +  L 


U2) 


PWAKe'7'  = 


x'  =  x.  +  L 

^Hr^   I    Q(x»,y)dx« 
a  u  J 

X1  =  X. 


U*) 


where  the  geometry  is  sketched  in  Figure  14.. 


33 


yj+1 


*\ 


HEAT  ADDITION 
ZONE 


WAKE 


X| 


Xi 


Xj+L 


-►X 


FIGURE  14.  GEOMETRY  OF  WAKE  DOWNSTREAM  OF  A  HEAT  ADOmON  ZONE. 
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For  a  point  inside  the  heat  release  region 


pWAKE(xI,y)  =  "  ^T^    I   Q(x',y)dx»  (U) 

a  u 


The  symbol  xT  is  defined  in  Figure  14. 
Heat  Addition  at  a  Plane;  A  Heat  Front 

Consider  a  distribution  of  line  sources  which  form  a  plane  at  angle  9 
to  the  flow,  as  illustrated  by  Figure  15.   Terminology  similar  to  that  used 
in  supersonic  wing  theory,  i.e.,  subsonic  and  supersonic  "leading  edge,"  is 
employed  to  identify  two  cases,  9  ^  p.   or  M  ^  1.   The  normal  Mach  number  M 
is  M  sin  9.   Since  sin  p.  =  l/M,  the  normal  Mach  number  is 

M  =  ^  (45) 

n   sin  u  ' 

You  can  see  from  an  examination  of  Figure  15  that  the  heat  source  can  influence 

flow  upstream  if  M  <T  1.   However,  if  M  >  1,  the  heat  source  does  NOT  influence 

flow  upstream  contrary  to  statements  made  by  Tsien  and  Beilock.      Also  it 

is  apparent  that  point  P  is  influenced  only  by  right-running  characteristics  for 

M  <  1,  whereas  P  is  influenced  by  both  families  of  characteristics  for  M  >  1. 

n    '  *  n 

The  perturbation  velocities  for  a  heat  front  will  now  be  obtained. 
Using  Equation  (22)  and  integrating  along  the  heat  front,  the  formula  for  u'  is 


ui(x  v)  =  _  (r  -  i)q 

U  KX'yj  2Ypp 


r 
5  Kx  -  s  cos  9)  -  p(y  -  s  sin  0)1  ds        (46) 


-oo 
where  s  is  distance  along  the  heat  front.   Consequently, 
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,  ,  (r  -  i)g 


-  1 


U'    = 


2yp£       f  sin  ©  -  coa  © 


J 


oo 
6 

-60 


[•" 


*  -  By 

cos  ©  -  £  sin  © 


ds  = 


{P) 


,    <Y  -  l?q - 

2YPp(f   sin  ©  -  cos  ©) 


Since   £  =    -yM     -  1,    it  is  also  true  that  £  =  l/tan  ju.      Using  this  fact 

ui  =    (T  -  l)q  tan  n  sin  u 
2yp  sin(ji  -  0) 


as) 


For  left-running  waves,  p.   is  positive;  and  for  right-running  waves,  p.   is 
negative.  For  a  subsonic  "leading  edge,11  ©  is  less  than  u.   Downstream  of  the 
heat  source  for  M  <  1,  Equation  (4-8)  indicates  u'  will  be  negative;  upstream 


of  the  source  u1  will  be  positive.   In  addition 
u        r 


u'  > 
ul 


u 


(for  M  <  1) 
n 


where  subscripts  u  and  d  indicate  upstream  and  downstream  respectively. 

For  a  supersonic  "leading  edge,"  which  is  illustrated  on  the  right-hand 
side  of  Figure  15,  the  left  running  Mach  wave  gives  a  uJ  which  is  negative. 
The  right-running  Mach  wave  gives  a  ui  which  is  also  negative  but  which  is  less 
in  absolute  value  than  uJ.  Subscripts  L  and  R  denote  left-  and  right-running 
waves.   Consequently, 

u-  =  u£  ♦  u£  <0 

Now  let's  look  at  v';  from  Equation  (21)  one  obtains  for  this  case: 

00 

6  |(x  -  s  cos  ©)  -  £(y  -  s  sin  0)1  ds  (49) 


r 


i(x  y)  =  (r  -  i)q 

K    f7J  2Yp 


J 


-00 
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Integration  of  the  6 -function  yields 

v,  =  (r  -  l)q  =  (Y  -  1)q  sin  p  ,      . 

2yp    P  sin  9  -  cos  9     2yp    sin  (9  -  ju)  wu; 

For  the  case  M  <  1,  v'  is  negative,  and  v'  is  also  negative. 


Furthermore 


v'  >  VJ 

u  \  '    \   d 


(for  M  <  1) 


Fo 


r  the  case  M  >  1,  v'  is  positive  while  v'  is  negative.   In  addition 
n       1j  xt 


V*  =  v. 


'  -  v'   >0  (for  M  >  1) 


v  =  -  e 
r      x 


d    L   1  R I  ^    N     n 

which  means  the  flow  is  turned  toward  the  line  source  much  the  same  as  an 

oblique  shock  wave.   For  M  <  1,  the  flow  turns  away  from  the  heat  front. 

A  natural  question  is,  "Is  there  any  special  direction  to  the 
perturbation  velocity  vector?'1  Let  v  be  the  resultant  perturbation  velocity 
vector.   It  is  given  by 

J  tan  u  sin  u   tan  u  sin  u\  .  ~T  A     sin  u       sin  u  \         /__  N 
cG(sin(9  -u)  *  sin(9  !  u)  j  +   ey>in(9  -   u)  '   sin (9  Tfi)        (5l) 

The  algebraic  signs  are  correct  for  the  case  M  >  1  which  means  9  >  u.   The 

symbol  C  equals  (y  -  l)q/2yp.   Let  a  be  the  angle  of  the  perturbation 

velocity  vector  relative  to  freestream  direction.   Then 

sin  u        sin  ju 

tan  a  =   fin<9  "> *in(Q  ^3  (52) 

tan  u  sin  u   tan  ja  sin  ju 

sin(9  -  }x)         sin(9  +  ju) 
Manipulation  of  Equation  (52)  leads  to 

1   sin  u  cos  9       1  /_-\ 

tan  o  =  -  7 c  — : — 5  =  -  t 5  (53; 

tan  fji   cos  u  sin  9     tan  9 

The  conclusion  from  Equation  (53)  is  that  the  velocity  vector  is  normal  to 

the  planar  heat  source.   Since  this  is  so,  one  can  transform  the  supersonic 

case  of  heat  addition  to  an  equivalent  ID  problem.   To  do  this,  one  uses 
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(54) 


u  u 

where  u  =  U  sin  0.   Proceeding  by  substitution  of  components  from 
Equation  (51)  into  Equation  (54.)  gives 


du        1     h. 

1  (55) 


U  sinQ    1  -  M2  c  T 
n  p 

which  is  precisely  the  influence  coefficient  result  when  properly  interpreted. 
See  Table  IV.   For  a  subsonic  "leading  edge,"  the  velocity  increases.   This 
means  that  the  flow  is  turned  away  from  the  planar  heat  source  for  M  <  1 
(analogous  to  Prandtl  Meyer  expansion)  and  toward  the  planar  heat  source  for 
M  >  1  (analogous  to  oblique  shock  wave). 

As  has  been  demonstrated,  one  of  the  ways  of  adding  heat  which  lends 
itself  to  analysis  is  by  means  of  a  heat  front.  The  example  chosen  by 
Tsein  and  Beilock     was  a  heat  front.   Oswatitsch     considers  other 
aspects  of  heat  fronts. 

Heat  addition  at  a  front  is  comparable  in  complexity  to  linearized  heat 
addition.   Combustion  is  a  volume  phenomenon.   However,  if  the  length  of  the 

heat  release  region  is  much  less  than  a  body  length,  a  front  is  a  reasonable 

(12) 
approximation.   Detonation  is  a  heat  front  of  a  very  special  nature. 
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A  heat  front  with  a  subsonic  leading  edge  located  next  to  a  wall  will  be 
preceded  by  an  oblique  shock.  The  oblique  shock  turns  the  flow  away  from  the  wall, 
and  the  heat  front  turns  it  back  parallel  to  the  wall.  Recall  that  for  a  subsonic 
leading  edge  the  heat  front  can  influence  the  flow  upstream.   The  corresponding  changes 
in  pressure  are  shown.   The  net  change  in  pressure  is  positive  although  the  heat 
front  decreases  pressure. 

A  heat  front  with  a  supersonic  leading  edge  will  be  followed  by  an  expansion 
fan  if  the  front  :'.s  adjacent  to  a  straight  wall.   Once  again  the  net  pressure 
change  is  positive.   The  heat  front  turns  the  flow  away  from  the  wall,  and  the 
fan  redirects  it  so  as  to  be  parallel  to  the  wall. 

Adding  heat  to  a  supersonic  stream  decreases  the  Mach  number.  When  sufficient 
heat  has  been  added,  the  normal  component  of  the  Mach  number  downstream  of  the 
heat  front,  M  „,  will  become  unity.  The  amount  of  heat  required  to  make  M  ?  =  1  is 
known  as  "critical"  amount  of  heat.  With  a  heat  front  having  a  supersonic  leading 
edge,  heat  can  be  added  until  M  _  =  1.   If  more  heat  is  added,  the  flow  illustrated  in 
Figurel6(b)  changes  to  that  shown  in  Figure  16(a).   The  oblique  shock  decreases 
the  flow  Mach  number  changing  the  leading  edge  of  the  heat  front  from  supersonic  to 
subsonic. 

Adding  heat  to  a  subsonic  flow  or  a  heat  front  with  a  subsonic  leading  edge 
modifies  the  flow  upstream.   One  cannot  specify  the  flow  immediately  ahead  of  the 
front.   Critical  heat  addition  does  not  have  the  same  impact  in  the  subsonic  case 
since  the  flow  upstream  will  be  modified  to  adjust  to  the  large  heat  addition. 

A  detonation  wave  can  be  thought  of  as  a  shock  wave  followed  by  a  combustion  zone 
or  heat  front.   In  Figure.  16(a)  if  the  oblique  shock  and  the  heat  front  had  the  same 
angle,  it  would  be  an  oblique  detonation  wave  provided  the  Mach  number  and  amount  of 
heat  released  were  correct. 
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FIGURE  16.  HEAT  FRONTS  AND  COMBUSTION  FANS. 
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In  the  Prandtl -Meyer  expansion,  conditions  are  constant  along  a 
characteristic.   Centered  combustion  fans  are  possible  with  conditions 
constant  along  a  ray.   However,  the  ray  is  not  a  Mach  line,  and  the 
normal  component  of  Mach  number  need  not  be  unity.  Combustion  fans  could 
occur  in  subsonic  as  well  as  supersonic  flow. 

If  M  ,  the  normal  component  of  Mach  number,  exceeds  unity,  the  fan  will 
turn  the  flow  away  from  the  wall  like  a  compression.   This  is  illustrated  in 
Figure  16(c)  on  the  left-hand  side.   If  M  is  less  than  unity,  the  flow  would 
be  turned  downward  as  shown  in  Figure  16(c)  on  the  right-hand  side. 

Starting  with  the  continuity  equation  in  cylindrical  coordinates  (r,f,z) 
and  setting  9(  )/8z  =  0  and  8(  )/9r  =  0,  one  can  derive  the  continuity 
equation  for  a  combustion  fan.   It  is 

dw?     d-tnp 

w  ♦  — -  ♦  w~ =0  (56) 

1   dp    2  dp  VPD; 

The  velocity  components  are  illustrated  in  Figure  17,   Heat  addition  does  not 

cause  an  irrotational  flow  to  become  rotational.   One  additional  equation  is 

the  condition  for  irrotational  flow  which  is 

dw, 

dp  -  "2  -  °  <^) 

The  rate  of  heat  release  controls  d-unp/df.   One  simple  combustion  fan  occurs 
for  dvnp/dg  =  constant.   For  this  case,  Equations  (56)  and  (57)  can  be 
readily  integrated. 

There  is  a  particularly  interesting  expression  that  can  be  derived  for 
the  pressure  coefficient  associated  with  a  heat  front  in  hypersonic  flow. 
Consider  a  heat  front  at  angle  0  relative  to  the  main  stream.   The  velocity 
components  and  angles  are  illustrated  in  Figure  18.   The  momentum  equation  is 


U2 


REFERENCE  DIRECTION 


FIGURE  17.  GEOMETRY  FOR  COMBUSTION  FAN. 
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FIGURE  18.  DEFINITION  OF  ANGLES  AND  VELOCITY  COMPONENTS. 


44 


Pi  +  pl^  =  P2  *  P2U2  (58> 


which  can  be  rearranged,  since  from  continuity  p,u.  ■  P2U? 


P2  "  Pl  =  plUl  -  P2U2  =  Pl^S  ~  V  (59) 


From  the  geometry  illustrated   in  Figure  18 

W2  *1   "  U2 


sinj  p  -  ©I  sin  6 


(60) 


Noting  that  11.  =  w,  sin  0  and  combining  Equations  (59)  and  (60)  leads  to 

P2  -  p,  =  Piwiw2  tan  ®  sin^  (61) 

If  M  =  H  sin  9  >>   1,  heat  addition  has  little  influence  on  w  or  p.   Hence 
n 

w,  <-  w_.   For  hypersonic  flow  Equation  (61)  becomes 


P2  -  Pn  =  2  2Piwi)tan  9  sin 6 


and  the  pressure  coefficient  is 


c  =  2  tan  8  sin6  (62) 

P 

This  can  be  compared  with  the  Newtonian  expression  which  is 

c  =  2  sin26  (61) 

P 

Equations  (62)  and  (63)  provide  sufficient  tools  to  explore  the  external 
burning  ramjet  shown  in  Figure  19.   Angle  5   is  shown  in  the  top  figure  and 
applies  to  Equation  (63). 
Heat  Addition  and  Vorticity 

In  this  section  it  is  shown  that  heat  addition  does  not  change  an 
initially  irrotational  flow  to  a  rotational  flow.   The  curl  of  the  velocity 
vector  should  be  zero.   For  subsonic  flow,  one  uses  the  equations  from  Tsien 
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and  Beilock    ;  and  for  supersonic  flow,  one  uses  Equations  (22)  and  (23). 

curl  7=  rjf^'  -f^j  (6A) 

Substituting  Equations  22  and  23  into  Equation  (64.)  yields 

curl  7=  7zl^±l2[-   2n(x  -  for)  (x  -  py)]  ♦  ^^[2nP(x  -  py)  (x  -  py)]|  (65) 

Recall  that  p  =  M  -  1.   In  the  above  analysis,  5   has  been  represented  by 


*(*  -  py)  =  nl«^  exp[-  n(x  -  Py) 


2 


(66) 


Since  curl  v  vanishes,  one  concludes  that  heat  addition  does  not  generate 

vorticity. 

Mass  Balance  with  Heat  Addition 

To  gain  additional  insight  to  diabatic  flows,  it  is  worthwhile  to  examine 
some  simple  flows  using  familiar  methods.   One  such  case  is  a  mass  balance. 
Consider  a  heat  addition  region  as  shown  in  Figure  20  with  a  control  volume. 
Half  of  the  control  volume  can  be  used  for  the  mass  balance.   This  half  has 
width  w.   Writing  a  continuity  equation  for  mass  flux  across  control  volume 
yields 

P*u«,w=  (P  +  P'V(L  -  pi)  *  2(p  ♦  p'K'pl  ♦  PoauM(w  -  1/2)   +  (Pco+  p^)u  1/2    (67) 

where  p'  and  v1  are  averages  of  perturbation  quantities  along  pi  illustrated 

in  Figure  20.   The  density  perturbation  in  the  wake  is  p'.   Equation  (67) 

w 

reduces  to 

pv'L  =  -  9^1/2  (68) 

when  second  order  terms  (e.g.,  p'v')  are  neglected.   This  can  be  rewritten  as 
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1       4 
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v»     p«  I 

—  =  -—  oT  (69) 

Substitution  of  Equations  (4.0)  and  (4.2)  in  (69)  gives  an  equality,  as  one 
expects. 

EXTERNAL  BURNING;  THRUST  MINUS  DRAG 
Drag — Subsonic  and  Supersonic 

In  subsonic  flow  of  a  perfect  fluid  around  a  body,  there  is  no  force  on 
the  body,  either  lift  or  drag.   Motion  of  a  body  through  a  fluid  would  require 
work  if  there  were  a  drag.   This  energy  would  appear  in  the  fluid.   For  a 
perfect  fluid,  the  perturbations  in  the  fluid  decay  rapidly  at  large 
distances  from  the  body.   If  the  body  has  been  moving  for  a  long  time  and  if 
there  were  drag,  there  should  be  significant  motion  of  the  fluid  at  large 
distances. 

In  supersonic  flow,  there  is  significant  motion  at  large  distances  even  in 
a  perfect  fluid.   The  waves  caused  by  body  motion  move  out  from  the  body  eventually 
decaying  to  sound  waves.   Far  from  the  body,  the  flow  consists  of  outgoing  sound 
waves.   These  waves  carry  away  energy.   Work  overcoming  drag  supplies  the 
energy. 

Relation  Between  Flux  of  Entropy  and  Drag 

(13) 
It  is  well  known  (See,  for  example,  K.  Oswatitsch.     )  that  the  flux  of 

entropy  through  a  control  surface  can  be  related  to  drag  of  a  body  within  the 

control  surface.   To  get  a  feeling  for  the  extent  of  the  waves  causing  drag,  a 

sample  problem  was  worked.   The  results  are  shown  in  Figure  21. 

Consider  a  20  wedge  at  Mach  2,  as  shown  in  Figure  21.   The  wave  geometry 

is  shown  and  has  been  calculated  using  the  method  of  finite  waves.   Width  of 
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POINT      JL/b  PERCENTAGE  DRAG 
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179. 

FIGURE  21.    FLOW    FIELD  NEAR  A  WEDGE  AT  MACH  2 
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the  base  of  the  wedge  is  b.   The  lateral  distance  to  a  point  on  the  shock 
wave  is  I.      Distance  t~   is  illustrated.   The  flux  of  entropy  in  and  out  of 
the  control  volume  (neglecting  the  wake)  is  equal  to  entropy  jump  across 
the  shock  wave  times  the  mass  flux.   The  first  expansion  wave  intersects  the 
shock  wave  at  point  2.   The  wave  between  points  1  and  2  accounts  for  25.2 
per  cent  of  wedge  drag.   The  shock  wave  between  points  1  and  3  accounts  for 
39.3  per  cent  of  the  wedge  pressure  drag.   It  is  necessary  to  go  out  179 
base  widths  to  obtain  99  per  cent  of  the  drag. 

The  conclusion  from  this  little  study  is  that  the  waves  are  important  to 
great  distances  from  the  body.   In  axisymmetric  flow,  the  bow  wave  would  decay 
more  rapidly. 
Linearized  Drag  Due  to  Waves 

A  wave  is  a  mechanism  whereby  energy  can  be  propagated.   The  energy  takes 
two  forms  in  the  linearized  or  acoustic  case.   These  are  flow  work  and  kinetic 
energy.   It  can  be  shown     that  an  integral  over  space  of  the  flow  work  equals 
the  integral  over  space  of  kinetic  energy.   As  an  example,  consider  a  half- 
diamond  airfoil  shown  in  Figure  22.   The  airfoil  is  moved  exactly  one  chord 
length  from  position  ABC  to  position  DEA,  and  the  undisturbed  fluid  is  at  rest. 
This  motion  does  work  Dc  where  D  is  the  drag  and  c  is  the  chord.   The  work  appears 
as  energy  within  the  waves.  Changes  in  the  fluid  motion  occur  normal  to  the 
waves.   The  disturbance  originating  at  A  now  appears  at  A1,  B  at  B',  C  at  C1,  etc. 
The  motion  within  the  wave  below  A'B'C  is  due  to  earlier  airfoil  motion  to  the 
right  of  position  ABC.   The  work  done  by  drag  moving  one  chord  length  appears 
within  the  volume,  V,  bounded  by  A'B'C'DKA.   Consequently 
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FIGURE  22.   DRAG  EVALUATION   USING  WAVE  ENERGY. 
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Dc  = 


[p  (U  ♦  u)2  +  v2" 


dV 


(70) 


where  u  and  v  are  perturbation  velocities.   When  there  is  heat  addition,  as 
illustrated  in  the  lower  part  of  Figure  22,  Equation  (70)  changes.   The  force 
on  the  body  is  now  thrust  minus  drag.   Heat  has  been  added.   There  is  a 
left-running  wave,  AG,  which  moves  into  the  fluid  above  the  airfoil.   Energy- 
is  added  to  the  fluid  bounded  by  ACG.   Heat  is  added  in  a  volume  bounded  by  HFC. 
Downstream  of  this  region  there  is  a  wake.   Note  that  in  the  region  EAC'B1  the 
mass  sinks  cause  u1  and  v1  to  be  positive*  whereas  the  heat  sources  cause  u" 
and  v"  to  be  negative.   In  region  EAC'B'  the  mass  sinks  decrease  p  whereas  heat 
increases  p.   The  energy  equation  associated  with  the  waves  now  is 


(T  -  D)c  +  (heat  added)  =  |    p  (U  +  u'  +  u")2  +  (y«  +  v")2 


dV 


(TL) 


pd(i)  * 


( energy )dV  ♦    ( energy )dV 
"V         WAKE  ~AGC 

Equation  (71)  is  more  conceptual  than  quantitative.   There  are  some  questions 
yet  to  be  answered.   Does  the  integral  of  flow  work  equal  the  integral  of 
kinetic  energy  for  waves  caused  by  heat  addition?   Incorporating  heat  addition 
into  the  waves  increases  complexity.   Can  this  method  of  drag  and  thrust 
evaluation  lead  to  useful  results  and  new  insight?  These  questions  will  be 


*We  use  the  notation  that  mass  sources  cause  perturbation  velocities  u'  and  v1 ; 
heat  sources  cause  perturbation  velocities  u"  and  v". 
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considered  in  the  section  Wave  Energy  Due  to  Heat  Addition  which  follows 

shortly. 

Linearized  Momentum  Relations 

In  two  dimensions,  x  and  z,  the  drag  formula  in  three  dimensions  is 
simplified.   It  becomes 


|  (p  -  Pfl9)dz  -   I  p(U  +  u)[(U  +  u)dz  -  wdx]  (72) 


D  =  - 

"c 

The  corresponding  lift  equation  is 


L  =    I  (p  -  Pa>)dx  +   i  pwfwdx  -  (U  ♦  u)dz]  (73) 

G  C 

where  C  denotes  a  closed  contour  surrounding  the  vehicle  or  wing.  For 
small  perturbations  due  to  mass  sources  and  sinks,  the  above  equations 
become 


D  =  pm     I  uwdx  +  2Pa    [h  +  w  Jds 


(74) 


2       2 
Note  that  drag  is  due  to  second  order  terms,  i.e.,  uw,  u  ,  and  w  .   First 

order  terms  drop  out.   For  lift  the  small  perturbation  form  is 


L  =  -  P«u    (udx  +  wdz)  =  P«ur  (75) 


■"J 


G 

In  a  following  section  the  equivalent  formulas  for  L  and  D  will  be  derived 
for  combined  mass  and  heat  addition.   First,  consider  drag  or  thrust  on  a  heat 
addition  zone. 

5U 


To  illustrate  the  application  of  the  heat  addition  equations,  let  us 
calculate  the  thrust  of  a  line  source  of  length,  I,   oriented  parallel  and 

normal  to  the  stream.   Equation  (n6)   is  taken  from  Chapter  7  of  Liepmann 

(15) 
and  Roshko; 


D  =  - 


(pu^iy^dA  -  1  pi^dA 
A 


0>6) 


In  their  notation  x,  =  x  and  x2  =  y.   The  control  volume  is  shown  in  Figure  23. 
If  Equation  (76)  is  negative,  one  gets  a  thrust.   Using  the  numbers  for 
surfaces  in  control  volume: 


D  =  - 


(pU)(-  U)dA  - 


(P.- 


CO 


W dA 


(U  -  u')(p_  -  P'WAVE)(-  1)(-  V)dA  - 


00 


(77) 


r 


J 


(U  -u')(pao-p'WAVE)(*D(v')dA 


The  integral  over  A.  can  be  integrated.   The  integral  over  A?  can  be  split 
into  two  parts,  and  the  part  with  p^  can  be  integrated.   Equations  (22)  through 
(25)  were  derived  dropping  terms  in  second  order  in  v'/U,  etc.   The  integrands 
of  A^.  and  A.  reduce  to  P^v'.   Consequently,  Equation  (77)  becomes: 
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D  =    -  p^U2^   ♦  Po/a3   ♦  (P'WAKS^^  "  P-Uv'd^ 


(78) 


We  now  need  expressions  for  p  '»,.«•«  and  v1,   which  we  obtain  from 

Equations   (25)  and   (23)   respectively.      Proceeding  with  p' 

I 

x,y)  (t  ;l)q6(y-y')l(x  -x')ds 

a  U 
0 


0>9) 


where  s  is  distance  along  the  source;  x  and  y,  the  coordinates  where  you  want 
to  evaluate  pA.vttJ  anc*  x1  and  y1,  the  coordinates  of  an  element,  ds,  of  heat 
source.   When  x>  x',  1=1;  since  we  are  evaluating  A~-integral,  x  will  be 
greater  than  x'.   Also  ds  =  dx1  for  the  case  at  hand. 

.i 


a  U 


dx»  =  (t  r^sfr-y) 

a  U 


(80) 


Substituting  Equation   (80)   into  A~-integral  of  Equation   (78),   we  get 


p  w2"  ■        ^^T^Ky  -  y)U2dy  .  i^^is 


WAKE 


a  U 


(81) 


"3  *3 

In  a  similar  manner,  an  equation  can  be  obtained  for  v1  from  Equation  (23): 

rl 


v'(x,y)  = 


(Vp1)q6[(x"x,) "  *<y-y,>]ds       (82) 


0 


Substituting  Equation  (82)  in  A^-integral  of  Equation  (^8)  yields 
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Pa>u(Y2rpl)q5[(X  "  X,)    "   f(y  "  y^]dx,d3C  (83) 


The  coordinate  system  is  chosen  so  that  the  heat  source  starts  at  x1  =  0  and 
ends  at  x'  =  t.   Also  y1  =  0.   The  coordinate  of  point  c  in  Figure  23  is 
(£yc>yc)>  anc*  ^e  coordinate  of  point  d  is  (Pyv  +  "Wy^)*  Consequently,  the 
limits  of  integration  are  as  shown 


]]{:L^rh      J         j  6[(x  -  x')  -  py5]dx'dx 


A  change  of  variable  helps  evaluate  the  integral;   let  z  =  x  -  py^.. 


roo       2yp 


6(2  -  x')dzdx'  =  PwU(r  ;ftp)q  (85) 

J       « 

0        0 

The  integral  in  Equation  (19)  is  in  the  standard  form  as  shown  in  the 
Appendix  so  that  it  is  readily  apparent  that  it  has  value  I.      Now  combine 
Equations  (78),  (81),  and  (85)  to  give 

D  =  (t  -  i)q*u  _  J^^iM]  (86) 

a        \    2yp    / 

2     , 
which  vanishes  since  a  =  YP/P»   Supersonic  heat  addition  in  an  unconstrained 

(no  walls)  flow  gives  zero  thrust  or  drag. 

Now  look  at  the  case  where  the  heat  source  is  of  length  t   but  normal  to 

the  streamlines.   This  is  illustrated  in  Figure  24..   Equations  (78),  (79),  and 

(82)  apply;  however,  ds  is  different  as  is  the  limit  of  integration.   Put 

the  origin  of  coordinates  at  the  center  of  the  heat  source.   Equation  (79)  becomes 
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in 

-"»  [  s 

U    J 

-1/2 


p       (  )  =  icz 


The  A~-integral  in  Equation  (78)  becomes 

.1/2 


8  (y  -  y')dy« 


(87) 


I"2^  I 


6(y-y')dy'dy  =  (^  =  ^*U 


(88) 


a 


A?  -  1/2 

since  the  integrations  over  y'  and  y  overlap  over  the  region  -  1/2   to  1/2, 
For  the  perturbation  velocity,  the  equation  for  v'  becomes 

1/2 


r 


v'(x,y)  = 


^^^Jasfx-p^-y'jjdy' 


(89) 


J-l/2 


where  y^  is  the  coordinate  for  plane  A,.   The  A^-integral  is 


r 


1/2 


PooUv'dy  =     I    (T  ""^pqPmUS^  "  P(y5  "  y')]dy'dx 


(90) 


A5  J-  1/2 


The  limits  of  integration  for  x  are  pCy^  -  1/2)   to  p(yc  *  1/2).   To  obtain  a 
standard  form  for  the  integrand,  we  must  eliminate  f  in  the  6 -function.   This 
we  do  by  letting  z  =  P(yr  -  y' )  and  dz  =  -  fdy'.   Equation  (90)  is  now 


*»   2yp£ 


8  (x  -  z)dzdx  =  p  U^r  ~  yq 


2ypP 


P(y5  +  *)  -  P(y5  -|) 


(91) 
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The  limits  of  integration  for  both  z  and  x  in  Equation  (91)  are  fCyv  -  t/2) 
to  P(y"c  +  V2).   Combining  Equations  (78),  (88),  and  (91),  one  finds 


D  =  2 "  2I I  =  °  (92) 


(Y  -  DqlU    /(y  -  Dqlpjh 

"2  

a        y    2yp 

which  is  what  we  expected. 

It  is  interesting  to  draw  the  streamlines  in  the  vicinity  of  the  heat 

sources;  see  Figure  25.   For  the  parallel  heat  source  there  is  a  conceptual 

difficulty  which  has  a  resolution.   The  "stagnation"  streamline*  C  goes 

straight  back  since  v1  is  an  odd  function  and  must  be  zero  at  y  =  0.   One 

can  think  of  a  streamtube  bounded  by  streamlines  B  and  D.  Continuing,  one 

can  shrink  this  streamtube  so  that  there  is  very  little  mass  flow  through  it. 

Yet  the  streamtube  diverges  at  angle  0.   A  finite  amount  of  heat  is  added  to  an 

infinitesimal  stream  yielding  a  finite  angle  0.  Strange.'  One  avoids  this 

difficulty  by  using  volume  heat  addition.   The  case  of  the  heat  source  normal 

to  the  flow  is  more  easily  digested.   Streamlines  F  and  J  are  acted  upon  by 

characteristics  of  a  single  family.   Streamline  G  passes  through  the  heat 

source  into  a  region  where  it  is  influenced  by  both  left-running  (gives  +  v') 

and  right-running  (gives  -  v')  characteristics.   In  this  region  the  streamline 

is  not  deviated  from  its  original  direction.   Streamline  G  then  moves  into  a 

region  where  only  left-running  characteristics  act  on  it.   It  is  deflected  by  an 

amount  0.   Finally  G  moves  into  the  wake  where  it  accelerates  and  assumes  its 

original  direction.   The  speed  along  G  is  (U,  U  -  2u',  U  -  u1,  U)  as  it  moves 

from  region  to  region.   Streamline  H  is  in  the  plane  of  symmetry;  as  a  result, 

it  does  not  deviate  from  freestream  direction. 


*A  body  will  stagnate  the  flow;  however  a  heat  source  does  not  do  so. 
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FIGURE  25.  STREAMLINES  IN  VICINITY  OF  PARALLEL 
AND  NORMAL  LINE  HEAT  SOURCES. 
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Wave  Energy  Due  to  Heat  Addition 

The  model  for  determining  wave  energy  due  to  heat  addition  is  shown 
in  Figure  26.   A  heat  addition  zone  BFGH  is  translated  its  own  Length  L  in 
time  L/ug,.   Its  new  position  is  AFEL.   The  heat  released  is  (hll)  (L/uJ) . 
Since  there  is  no  force  on  the  heat  addition  region,  no  work  is  done  moving  it 
from  EFGH  to  AFEL.   The  question  of  forces  on  heat  addition  zones  was 
discussed  in  the  preceding  section.   All  of  the  heat  energy  released  must 
appear  in  the  waves  and  in  the  wake. 

In  the  time  L/ua,,  the  disturbance  at  J  has  moved  to  point  K.   The 
disturbance  at  F  has  moved  to  B.   The  dashed  lines  AC  and  BD  would  be  stream- 
lines if  EFGH  were  stationary  and  the  fluid  moved  by.   The  energy  content  of 
space  bounded  by  ANCEL  is  the  same  as  that  inside  FJ1-:HE.   Part  of  the  energy 
added  appears  in  the  region  ABDC.   This  is  the  wave  energy.   The  remainder 
appears  in  the  wake  in  the  region  EFGH.   From  B  to  P,  work  is  done  on  the 
fluid  to  compress  it.   From  K  to  D,  work  is  done  by  the  fluid  as  it  expands. 
The  flow  work  cancels.   We  are  now  ready  to  write  the  energy  equation. 

h-tL  ._  2  [kinetic  energy]    /  energy  in  \  /q,x 

u   "    \         in  ABDC    /    \   wake  EFGH  J  vv  ; 

The  factor  2  comes  from  the  fact  that  waves  go  both  up  and  down  in  Figure  26. 
The  kinetic  energy  in  ABDC  is 


2Pco  I      undxd^  "  |p- 


L  cos  u 


2 
un  sin  p.   djdn  (94-) 


J 


x  "y  0   0 

where  dx  and  dy  have  been  replaced  by  a  coordinate  system  x  =  £  and  y  =  i^  sin  ji. 
The  symbol  u  represents  the  perturbation  velocity  normal  to  the  wave.   It  is 
e  qual  to 

63 


B 


----- /M 


FIGURE  26.   WAVE  ENERGY  DUE  TO  HEAT  ADDITION. 
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u2  =  u'2  +  v'2  (95) 

n  ' 

where  u'  and  v'  are  given  by  Equations  (39)  and  (4-0)  respectively.   The 

magnitude  of  v'  has  been  plotted  in  the  upper  portion  of  Figure  20.   Inserting 

Equation  (95)  into  Equation  (94-)  and  recognizing  the  functional  form  of  v1 

(and  u' )  as  shown  in  Figure  20,  one  arrives  at 


2KE  in  ABDC 


,-a   j[(jf)  ■  <-* 


LPoosin  u  cos  u  I  \\  [ft]      '  \JZ 


r 


fdf  +  (u'2  ♦  v'2) 
mm 


dj  (96) 

pi 


The  value  of  v1   has   been  indicated  in  Figure   20.      It  is   the  maximum 
m 

perturbation  velocity.      Integrating  Equation  (96)  yields 

/      2KB  in  ABDC     \   .   2,     2  2)l  +  (     2  ,       2)(     .  ^  . 

/        3     m  m   /r  v  m  m  v   ' 


Lp  sin  a  cos  u 


(9*) 


=  (u'2  +  v'2)(L  --fr) 
m         m  3 


The  energy  in  the  wake   is 

(ZTw$  -  vP^Tw  -  v  ■  <vPv*£  ^ 

The  distance  FG  is  the  length  L.   According  to  small  perturbation  results, 

the  wake  has  freestream  velocity  and  pressure  but  reduced  density  and  increased 

temperature.   Knowing  these  facts  and  using  p  =  pRT,  Equation  (98)  becomes 


/ 


energy  m\ 


=  c 


TLlp«  (99) 


^wake  EFGHj    p  oo     *w 

The  density  perturbation  in  the  wake  p1  can  be  obtained  from  Equation  (44.). 

w 

Substituting  Equations  (39)  and  (4.0)  into  Equation  (97),  substituting 
Equation  (4-4.)  into  Equation  (99),  and  combining  the  results  in  Equation  (93) 
yields 
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QlL^       P<X)L  U 
=  (L  -  — ) 

9  * 


(y  -  l)Qt 


2yp 


2 


2i 


Y   -  1 


♦  c  T  L^Q-r-  (100) 

p  »  4 

r  oil 


WAVE  ENERGY  WAKE  ENERGY 

From  Equation   (100 )    it  is  apparent  that  the   energy  appearing   in  the  waves   is 

2     2 

second  order.   The  wave  energy  was  based  on  u'   +  v1  which  is  second  order. 

mm 
2 

The  wake  energy  reduces  to  Q(-L  /u^,.   If  Equation  (100 )  is  to  be  an  equation, 
the  wave  energy  must  be  neglected. 

There  are  some  important  conclusions  from  this  exercise.   First,  the 
wave  energy  is  second  order.   Second,  most  of  the  energy  due  to  heat  addition 
appears  in  the  wake.   Energy  radiated  away  by  waves  is  small  compared  to 
energy  deposited  in  the  wake. 

If  one  examines  the  linearized  momentum  equations  in  two  dimensions,  as 
given  in  the  section,  Linearized  Momentum  Equations,  one  finds  drag  is  related 
to  an  integral  of  second  order  terms.   As  shown  in  the  discussion  of  Figure  21, 
drag  can  be  related  to  the  waves  generated  by  motion  of  a  body.   A  conclusion 
is  that  wave  energy  resulting  from  either  drag  or  heat  addition  is  second  order 
in  the  perturbation  quantities. 
Control  Volume  Approach  to  Thrust  Minus  Drag  with  Heat  Addition 

As  an  example  of  a  two-dimensional  planar  flow,  examine  Figure  27.   A 
half  diamond  has  a  region  of  volume  heat  addition  indicated  as  a  shaded  area. 
The  dotted  lines  outline  an  image  heat  addition  region  above  the  shaded  region. 
Part  of  a  control  volume  boundary  is  shown  as  a  dashed  line.   Mach  lines  are 
drawn  from  the  body  to  the  boundary  of  the  control  surface.   Below  the  boundary 
are  plotted  the  perturbation  quantities  for  that  location  on  the  boundary.   At 
the  very  bottom  are  the  pressure,  source  strength,  u1  and  v1.   These  are  constant 
to  mid  chord  and  then  switch  to  values  opposite  in  sign. 
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The  perturbation  quantities  due  to  heat  addition  are  plotted  about  the  right- 
hand  side  of  mid-page.   Heat  addition  increases  pressure,  slows  the  flow,  and 
deflects  it  downward.   If  the  two  plots  were  to  scale,  the  positive  pressure 
increment  due  to  heat  addition  would  exceed  the  negative  increment  due  to  flow 
expansion  (mass  sinks).   The  body  would  have  positive  thrust  minus  drag. 

The  magnitudes  of  the  perturbations  due  to  heat  addition  are  proportional  to 
the  length  of  the  characteristic  (Mach  line)  passing  through  the  heat  addition 
regions  both  real  and  image.   Note  that  p',  p",  u1,  u",  v',  and  v"  are  of 
opposite  signs. 

Let's  look  at  drag  and  lift  contributions  due  to  the  portion  of  control 
volume  shown.   In  Equation  (72)  dz  is  zero.   Hence  the  drag  contribution 

D  =  I  p(U+u'  +  u,,)(wt  +  w")dx  (101) 

Multiplying  the  terms  yields 


D  =    (pUw 


•  +  pUw"  +  pu'w1  +  pu'w"  ♦  pu"w'  *  pu"w")dx  (102) 

12      3      4      5      6 
The  following  summary  gives  the  contributions  (+),  0,  or  (-). 


Term           1 

2 

3 

U 

5 

6 

0  to  c/2      (-) 

0 

(♦) 

0 

0 

0 

c/2  to   c      (+) 

(-) 

(+) 

(-) 

(-) 

(♦) 

The  integral  of  pUw'  over  a  full  chord  gives  zero;  however,  the  integral  of 
pUw"  makes  a  large  thrust  contribution.   Terms  3  and  6  give  drag,  whereas  terms 
2,  L,   and  5  give  thrust. 
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Lift  is  obtained  from  Equation  (73).   Once  again  for  the  segment  of 
control  volume  illustrated,  dz  is  zero.   The  lift  contribution  becomes 

|  (p«  +  p")dx  +   I  p(w'2  +  2w'w'»  +  w»2)dx  (101) 

12       J         3      U  5 

The  following  summary  gives  the  contributions  of  the  various  terms. 


Term           1 

2 

3 

U 

5 

0  to   c/2      (+) 

0 

(+) 

0 

0 

c/2  to   c      (-) 

(+) 

(*) 

(-) 

(+) 

The  pressure  p1  averages  to  zero,  whereas  prt  gives  a  lift.   The  cross  term 
w'w"  reduces  lift. 

To  optimize  the  heat  addition,  the  airfoil  shape  must  be  expressed 
analytically.   Heat  addition  also  needs  to  be  described  by  some  function  of 
x  and  z.   Heat  release  need  not  necessarily  be  uniform. 
Thrust  from  Cycle  Point-of-View 

The  thrust  due  to  a  propulsion  device  is  given  by 

F  =  m(w  -  wj  *  Ae(pe  -  VoD)  (104) 

If  one  evaluates  the  equation  at  the  Treff tzplane,  the  pressure  term  may  or 
may  not  drop  out.   This  is  one  item  which  Oswatitsch     discusses.   In 
Equation  (104.),  F  is  thrust,  and  m  is  mass  flow  rate  of  fluid  which  receives 
a  velocity  increase  w  -  w^j  over  the  freestream  velocity  w^. 

Consider  three  flows  as  shown  in  Figure  28.   Heat  is  added  at  stagnation 
conditions  and  at  some  Mach  number  other  than  zero.   The  pressures  at  the 
exhausts  are  all  equal,  i.  e. ,  p  =  p.  =  p'  =  plj.   The  Mollier  Chart  which  depicts 
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w; 


w; 


1 


INFINITE  CROSS  SECTION 


Wa 


J 


dq 


W2 


(a)    ISENTROPIC  WITH  NO  HEAT 
ADOITION;THIS  IS  A  REFERENCE 
FLOW. 


(b)    ISENTROPIC  WITH  HEAT 
ADDITION  AT  STAGNATION 


w; 


W< 


3        4 


dq 


(c)  ISENTROPIC  TO  STATION  2,  ADD 
HEAT  2  T03;  ISENTROPIC  3  T04. 


FIGURE  28.  THREE  FLOWS  WITH   TWO  DIFFERENT  WAYS  OF  ADDING  HEAT. 
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the  three  flows  is  shown  in  Figure  29.   It  is  apparent  from  energy  considerations 
that 


dq  +  ¥2    =  hi   +  {\  ~   h2}  (105) 


Continuing 


|(v*  -  w^2)  =  wdw  =  dq  -  (h  -  hp  (106) 


Now 


dw  >  0 


if  dq  >  (h,  -  hj)  =  dh  (1CP) 


w,  >  w' 


The  slope  of  a  p  =  constant  curve  is 


r5)  =  T  <108) 


'p 

Combining  Equations  (106)  and  (108)  leads  to 


wdw  =  dq  -  dh  =  dq  -  ||  ds  (109) 

P 

The  change  in  entropy  of  the  flow  is  given  by 

ds  =  |a  (110) 

q 

A  subscript  q  is   used   to    indicate  this    is   the   temperature  at  which  heat   is 
added.      Equations    (109)   and    (110)   yield 

wdw  =   dq(l    -  |~)  (111) 

q 

If  one  compares  flow  (a)  with  flow  (b)  of  Figure  23,  then  the  temperature 

ratio  in  Equation  (111)  is  T'/T--,.   If  one  compares  flow  (a)  with  flow  (c),  the 

temperature  ratio  is  T'/T  .   To  have  dw  >  0,  one  must  have  T  <"  T  . 
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hT+dq... 
hT 


ENTROPY 


FIGURE  29.  MOLLIER  DIAGRAM  FOR  HEAT  ADDITION. 
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As  a  numerical  example,  consider  flight  at  M  =  0.9.   Heat  is  added  at 
M  =  0  in  flow  (b)  and  at  M~  =  0.5  in  flow  (c).   Enough  heat  is  added  to 
accelerate  the  flow  so  that  NL  =  0.6.   The  same  amount  of  heat  is  added  in 
flow  (b).   The  results  are  that  v./v£  =  1.055  and  w"/v'  =  1.091. 

Equation  (ill)  is  not  usually  applicable  to  supersonic  or  hypersonic 
flow.   The  entropy  change  in  the  exhaust  was  attributed  solely  to  heat 
addition.   In  supersonic  flow,  shock  waves  cause  entropy  changes.   Further- 
more heat  addition  may  alter  the  shock  wave  geometry. 
Thrust  and  Drag  in  Various  Flow  Regimes 

For  subsonic  flow,  the  pressure  term  in  Equation  (104 )  is  zero.   For 
supersonic  flow,  the  pressure  term  if  evaluated  in  the  Trefftzplane  is  not 
zero;  however,  it  is  extremely  small.   To  see  that  the  pressure  is  not  zero, 
examine  the  flow  illustrated  in  Figure  30.   The  entropy  in  the  Trefftzplane 
is  greater  than  ambient;  ds  >  0.   If  dp  were  zero,  then  pwdw  equals  -  Tds 
so  that  dw  <  0.   Entropy  is  related  to  p  and  T  by  Tds  =  c  dT  -  dp/p  =  c  dT. 
Consequently  dT  >  0.   From  p  =  pRT,  it  follows  dp/p  <  0.   The  continuity 
equation  pwA  must  be  satisfied.   As  shown  in  Figure  30,  the  streamtube  area  A 
does  not  change  significantly.   The  changes  in  both  p  and  w  cannot  be  negative 
and  still  satisfy  continuity.   If  dp  >  0,  then  dw  continues  to  be  less  than 
zero  and  dT  >  0.   Also  dp/p  >  dT/T.   Consequently  dp  >  0. 

To  account  for  the  drag,  the  area  A  must  be  quite  large.   (See  Figure  21. ) 
The  streamtube  influenced  by  a  propulsion  system  will  be  small  as  shown  in 
Figure  30.   It  is  an  area  of  approximately  A  in  Equation  (104.).   The  pressure 
at  the  Trefftzplane  is  due  mainly  to  drag  and  is  spread  over  area  A,  albeit  not 
uniformly.   The  pressure  term  in  the  thrust  equation  is  applicable  to  area  A  . 
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SUBSONIC 


STREAMTUBE 
FOR  PRGPULSIO 
OF  AREA  *>  Ae 


DRAG  IS  DUE  TO  DIP  IN 
VELOCITY 


PRESSURE   IN  WAKE 
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FOR  STREAMTUBE 
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}Ao 


FIGURE  30.    DRAG  IN  DIFFERENT  FLIGHT  REGIMES. 
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Since  thrust  equals  drag,  the  pressure  term  will  be  of  the  order  of  A  /A.   It 

e 

can  be  neglected  for  subsonic,  transonic,  and  supersonic  flow. 

For  slender  bodies  in  hypersonic  flow,  dw/w  is  small.   If  the  thickness 
ratio  is  T  =  ( thickness )/( body  length),  then  dw/w  is  of  the  order  t.      Heat 
addition  also  has  little  influence  on  dw/w,  a  fact  that  can  be  verified  by 

examining  Figure  7.3  in  Volume  I  of  A.  H.  Shapiro's  The  Dynamics  and  Thermo- 

(7) 

dynamics  of  Compressible  Fluid  Flow.     In  hypersonic  flow,  heat  addition 

causes  large  changes  in  pressure  and  static  temperature. 

In  hypersonic  flow,  the  shock  waves  are  close  to  the  body.   The  drag  is 
concentrated  in  a  small  area,  A-.;  see  Figure  30.   Thrust  producing  devices 
would  influence  an  area  comparable  to  A^.   Consequently,  the  pressure  term  in 
Equation  (104.)  must  be  retained  as  a  significant  factor. 

One  can  calculate  forces  on  a  body  by  suitable  integration  of  pressure. 

Pressure  in  excess  of  ambient  on  forward-facing  surfaces  causes  drag;  define 

this  as  pressure  drag,  D  .   Pressure  less  than  ambient  on  rearward -facing 

surfaces  causes  drag:  define  this  as  suction  drag,  D  .   The  ratio  D  /D  is 

to'  b*   s  s  p 

shown  in  Figure  31  as  a  function  of  Mach  number.   Suction  drag  is  very 
important  in  subsonic,  transonic,  and  supersonic  flow.   However,  as  M 
becomes  large  and  enters  the  hypersonic  region,  D  becomes  insignificant. 

9 

In  fact,  Newtonian  theory  assumes  ambient  pressure  p  is  zero  with  the  result 

D  =  0. 
s 

A  propulsion  device  will  increase  pressure  on  rearward  surfaces  as  shown 
in  Figure  31.   Heat  addition  must  increase  pressure  so  that  the  area  under  the 
curve  labelled  "with  heat"  equals  the  area  under  the  curve,  D  .   In  supersonic 

flow,  the  thrust  device  simultaneously  cancels  a  large  D  and  balances  D  . 

'  J  6   s  p 
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MACH  NUMBER 


SUPERSONIC 


DISTANCE   ALONG  BODY 


HYPERSONIC 


WITH   HEAT 


AMBIENT  p 


DISTANCE    ALONG  BODY 


FIGURE  31  CONTRIBUTIONS  TO  DRAG  AND  THRUST. 
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In  hypersonic  flow,  the  thrust  device  does  the  same  thing  except  that  in 
hypersonic  flow  D  is  negligible.   Heat  addition  in  hypersonic  flow  balances 
D  but  counteracts  only  an  insignificant  D  .   This  has  implications  con- 
cerning propulsion  efficiency. 

One  more  point  is  that  in  hypersonic  flow  D  is  the  major  drag  con- 
tribution, and  this  can  be  related  to  the  bow  shock  wave.   In  supersonic  flow, 
D,  is  important  also  so  that  the  waves  originating  at  the  rear  end  of  the  body 
or  at  the  rearward-facing  surfaces  are  of  equal  importance  to  the  bow  shock 
wave  in  calculating  drag. 
Contrast  in  Thrust  at  Supersonic  and  Hypersonic  Speeds 

Consider  a  body  immersed  in  a  flow  as  shown  in  Figure  30.   The  force  on 
the  body  in  the  direction  of  the  freestream  is 


=   I  [(pu2  -  Pau~)  *  (p  -  pjjdydz  (112) 


F  = 

A 

If  there  is  no  energy  addition  and  if  A  is  the  area  of  the  large  streamtube, 
F  equals  D.   (it  has  been  assumed  that  expansion  waves  have  made  pressure 
increments  negligible  along  the  streamline  from  a  to  b  in  Figure  30. )   If  there 
is  energy  addition  to  the  stream  and  if  A  is,  once  again,  the  Trefftzplane  in 
Figure  30,  Equation  (112)  yields  thrust  minus  drag.   If  one  evaluates 
Equation  (112)  using  the  small  streamtube  in  Figure  30  with  flow  m,  then  F  is 
approximately  thrust.* 


*It  is  not,  of  course,  possible  to  separate  out  the  thrust  accurately  using 
Equation  (112).   If  it  were  possible,  how  much  more  simple  the  propulsion 
integration  problem  would  be. 


77 


The  continuity  equation  is,  neglecting  fuel  added, 

0  =   \  (pu  -  p^ujdydz 
The  energy  equation  in  integral  form  is 


(in) 


o  = 


2    2 
u  *  v 


u 


pu(h  + 


00 


)  -  P0Ou-,(hfiO+—  *  q) 


dydz 


(114) 


In  Equation  (114.),  q  is  the  heat  added  per  unit  mass,  and  v  is  the  component 
of  velocity  normal  to  the  freestream.   By  combining  the  preceding  equations 
in  the  following  way 


u 


80  , 


(u^  x  thrust)    +•  (hM  +  q  -  -5-)   x  continuity  -  energy  =   u^F 
one  obtains 


v- 


pu(h.  -  h   +  q)    ♦  u   (p   -   p)    -  pu 


f  ^2   l    2 

(u    -   U„)        *  V 


so    *  '*> 


dydz 


(115) 


Substituting  h  =  e  +  p/p  into  Equation  (115)  yields 


usF  = 


ipu^q  -  (e  -  ej 


-  p(u  -  uj  -  pu- 


(u  -  uj   +  V 


dydz      (116) 


There  are  three  terms  in  the  integrand  of  Equation  (116).   The  term  pu(e  -  e^,) 
is  the  flux  of  energy  thrown  away  in  the  jet.   The  term  pu  jq  -  (e  -  e^J  is  the 
energy  converted  to  work  or  jet  kinetic  energy.   The  term  p(u  -  u^)  is  the  flow 
work,  and  the  last  term  in  the  integrand  is  the  flux  of  kinetic  energy  of  the 
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jet  relative  to  the  stationary  surroundings.   The  left-hand  side  is  the  rate 
at  which  the  thrust  does  work  in  overcoming  drag. 

For  transonic  and  supersonic  flow,  all  terms  are  important.   In  hypersonic 
flow,  (u  -  UpJ/u^,  v/u^jp  <^  1.   Changes  in  velocity  due  either  to  shock  waves  or 
heat  addition  are  small.   However  Ap/p  ,  AT/T^,  or  Ap/p^  are  of  the  order  of 
unit}'"  or  larger.   Based  on  these  comments,  one  can  neglect  the  jet  kinetic 

energy  in  hypersonic  flow.   In  hypersonic  flow,  the  term  p(u  -  ua)  varies  as 

2 

l/M  as  will  now  be  demonstrated.   Rewriting  the  term  as 


p(u  -O  -  P  "J1   +^)(—  "  D 

2 
one  can  introduce  the  Mach  number  dependence.   In  hypersonic  flow,  Ap/p  ~  K 

where  K  is  the  hypersonic  similarity  parameter  TM.   The  velocity  perturbation 

2 
varies  as  T  as  discussed  previously.   Substitution  of  these  quantities  into 

the  flow  work  term  yields 

P(u  -  uj  =  ^2(1  ♦  K2)(K2)  (117) 

Compared  to  the  heat  addition  term  in  Equation  (116),  the  term  p(u  -  u^)  is 
small  for  M  »  1  and  K  =  0(1 ).*  Equation  (116)  reduces  to 

F  =    p[q  -  (e  -  ejjdydz  (118) 

JA 
which  is  valid  for  K  »  1  and  K  =  0(l).   In  hypersonic  flow,  Equation  (118)  tells 
us  that  the  thrust  is  simply  due  to  the  heat  added  less  the  amount  of  internal 
energy  thrown  away  in  the  exhaust.   In  supersonic  flow,  it  is  more  complex;  see 
Equation  (116). 


•Requiring  K  =  0(1)  is  not  a  severe  restriction. 
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APPLICATIONS  OF  EXTERNAL  BURNING 

In  the  preceding  sections  some  of  the  analytical  tools  for  understanding 
external  burning  have  been  discussed.   It  is  worthwhile  to  look  at  some 
possible  applications  of  external  burning.   Four  applications  will  be  examined 
briefly  in  the  following  section;  these  are  forces  on  a  planar  airfoil  due  to 
heat  addition,  transonic  boattail  and  base  pressure  alleviation,  base  pressure 
modification  in  both  the  planar  and  axisymmetric  cases,  and  spin  recovery  of 
aircraft. 

Before  discussing  the  applications,  a  model  for  pressure  rise  due  to  heat 
addition  near  a  flat  plat  will  be  discussed.   This  mode]  uses  some  of  the 
results  from  EXTERNAL  BURNING;  BASIC  EQUATIONS  AND  SOLUTIONS  and  from 
EXTERNAL  BURNING;  THRUST  MINUS  DRAG.   The  model  is  a  refinement  of  a  similar 
model  by  Billig.  l6^ 
Model  for  Heat  Addition  Adjacent  to  a  Flat  Plate 

Figure  32  illustrates  the  geometry.   There  are  three  zones.   Zone  I  is 
the  freestream  which  is  supersonic.   Zone  II  is  the  heat  addition  region 
separated  from  Zone  I  by  a  heat  front.   Zone  II  is  bounded  by  the  wall,  heat 
front,  and  streamline  c.   Zone  III  is  downstream  of  the  oblique  shock  com- 
prising the  region  below  streamline  c.   The  heat  addition  turns  the  flow  by 
an  angle  8.   The  streamline  deflection  angles  and  static  pressures  in  regions 

II  and  LTI  must  be  equal  along  streamline  c.   In  the  following  development, 

2 
terms  in  0  are  dropped  as  being  negligible. 

The  deflection  angle  is  approximately: 

9  =  -2 i  (119) 
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The  continuity  equation  is 

PlVl  =  P3V3  =  p3U3(yl  *  9x)  =  P2V2  (l20) 

Across  the  oblique  shock  there  is  a  pressure  increase  equal  to 


P2  -  P-l      Y^©      YMJ© 


(121) 


pi      V^  "1   P 

Equation  (121 )  defines  £.   Define  Q1  as  the  heat  addition  per  unit  mass  of 
air  flowing  in  the  streamtube  defined  by  streamlines  a  and  c.   The  energy 
equation  is 

°pTTl  +  «fl  =  °pTT1  (122> 

or  as 

°pTT2  *  «&J  -    =pTT1  (122) 

(15) 

Using  Equation  (7.14)  from  Liepmann  and  Roshko,     one  can  derive  the 
following  momentum  equation 

"  P2V2  +  P3U3y3  "  P272  +  p3y3  "  P20X  =  °  (l23) 

Combining  equations  (119),  (120 ),  and  (123)  yields 

P2U2y2(u3  "  V  *  (p3  "  P2)(y2  *  ^  =  °  (12i^ 

A  solution  exists  for  u,  =  u~  and  p,  =  p?.   Let's  examine  the  consequences 
of  that  solution.   From  the  continuity  equation  and  equation  of  state 

P2U2y2   p3u3y3 


T2       I 


(125) 


Since  u^  ■  u^  and  p^  =  p   this  reduces  to 


T,   y, 

r "  T  (126> 
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The  energy  equation  becomes 


CpT2   *  2U2   +  Q23  =    %h   +  \A 
or  (127) 

c  T„    *  Q'     =   c  T, 
p  2  23  p  1 


since  u-  =  u  .   From  Equations  (126)  and  (l27) 

J2  p  2 

In  the  notation  of  the  influence  coefficients  of  Table  IV,  Equation  (128) 
can  be  rewritten  as 

dA   dh 

T=T?  (129> 

A    c  1 

P 
In  view  of  the  equality  expressed  by  Equation  (129),  the  influence 

coefficients  give  the  result  dp/p  =  0  and  du/u  =  0. 

The  heat  added  by  the  heat  front  must  increase  pressure  in  the 

streamtube  a-c  by  an  amount  equal  to  pressure  rise  across  the  oblique  shock. 

Using  this  fact 

dp   YM?9   YM?Q' 

p    f    Pc^ 

Combining  Equations    (119)  and    (128),    the  heat  added  in    Zone   II   is 

aj,  -  cpT2ex/yi  (in) 

The  difference  between  T^/T,  and  1  is  the  order  of  0.   Hence 

Q£3  =  c^QVy-L  (132) 
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Total  heat  added  to  turn  the  flow  by  an  angle  9  is 

«4l  "  «12  +  «23  =  9(1  *  £>  (111> 

The  heat  addition  Q',  turns  the  flow  by  an  angle  9  and  increases  pressure 

by  the  amount  given  by  Equation  (1^0). 

Forces  on  a  Planar  Airfoil  in  Supersonic  Fl ight  Due  to  Heat  Addition 

A  problem  which  has  been  studied  by  several  investigators  is  external 

(17) 
burning  near  a  two-dimensional  wing.   One  such  study  is  that  of  Mager. 

Mager  uses  the  linearized  heat  addition  formulas,  Equations  (22)  to  (25),  to 

find  the  forces  on  an  airfoil  in  supersonic  flight  with  heat  addition  in 

adjacent  streamtubes. 

Consider  an  area  bounded  by  x.  and  xf  in  the  x-direction  and  by  a 

surface  and  a  parallel  side  at  distance  h..   See  Figure  ?3.  Within  this  area, 

add  heat  Q.  From  Equation  (4.1)  the  pressure  at  the  surface  is 

p'  /y   -  1   QM  \     f  (Y  -  DQMh 

7 =  \—  IS]  J  sin  P  ds  =  -HF—  (1U) 

Where  does  the  factor  2  come  from?  There  must  be  an  image  heat  source  above  the 

surface;  otherwise,  there  would  be  a  v1  with  flow  through  the  wall.   Combining 

Equation  (134.)  and  the  definition  of  H  [h  is  defined  in  the  paragraph 

preceding  Equation  (3l).J  gives 

p^     (y  -  l)HhM 

p"  =  hi(xf  -  xi)afp  (1?5) 

In  Equation  (135)  the  quantity  h  varies  from  0  at  x  =  x.  to  h.  at  x  =  x.  +  h.  cot  u. 
It  remains  constant  to  x  =  x„,  and  then  h  decreases  linearly  to  zero  at 
x  =  x„  +  h.  cot  ji.   This  is  illustrated  in  Figure  34.. 


at 


SURFACE 


HEAT  H  UNIFORMLY 
RELEASED  IN  RECTANGLE 


FIGURE  33.    HEAT  ADDITION  ADJACENT  TO  A  SURFACE  IN  SUPERSONIC 
FLOW. 
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FIGURE  34.    PRESSURE  DISTRIBUTION  DUE  TO  HEAT  ADDITION 
ADJACENT  TO  A    SURFACE. 
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It  is  of  interest  to  examine  streamlines  near  the  surface  in  the  region 
of  heat  addition;  see  Figure  35.   The  streamline  passes  along  ABCD.   It  is 
shown  dashed  with  an  exaggerated  slope  below.   The  streamline  is  straight 
from  the  time  it  enters  the  heat  addition  region  until  it  arrives  at  A. 
This  is  so  since  aA  =  a 'A.   From  A  to  B  the  right-running  characteristic  grows 
in  length  relative  to  aA.   Hence  the  streamline  is  curved  from  A  to  B.   From 
B  to  C  the  difference  in  length  b'B  -  bB  is  constant;  hence  the  streamline  is 
straight  but  at  an  angle  to  the  mainstream.  From  C  to  D  less  and  less  of  the 
right-running  characteristic  passes  through  the  image  heat  addition  region, 
and  v1  becomes  smaller  and  smaller  until  finally  the  streamline  is  horizontal 
at  D.   From  C  to  D  the  streamline  is  obviously  curved. 

When  the  streamline  is  moved  closer  to  the  surface,  as  is  the  case  in 
the  lower  part  of  Figure  35,  the  deflection  is  considerably  less.   It  is 
obvious  that  the  streamline  on  the  surface  will  not  be  deflected  since 
aA  =  ~aJi  =  bB  =  cC  =  c^C  =  d*B  =  dD  for  this  case. 

Comparing  Figures  25  and  35 »  one  gains  some  insight  to  the  quandry 
discussed  in  connection  with  Figure  25.   When  the  heat  addition  is  confined  to 
a  mathematical  line,  as  in  the  top  of  Figure  25,  one  has  difficulty  because  of 
streamline  deflection  near  the  source.   However,  when  the  heat  is  distributed 
as  in  Figure  35 >  there  is  no  problem. 

Based  on  Figure  34.  and  Equation  (135),  it  is  apparent  that  one  has  in  hand 
the  tools  necessary  to  study  the  forces  on  bodies  due  to  heat  addition,  at 
least  in  the  linearized  case.   Refer  to  Mager's  paper  for  performance  of 
an  airfoil  with  heat  addition. 
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FIGURE  35.    STREAMLINES  IN  REGION  OF  HEAT  ADDITION. 
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Alleviation  of  Boattail  Drag  for  a  Plug  Nozzle 

In  the  Mach  number  range  0.8  ^  M„  ^1.2,  which  is  the  transonic  regime, 
the  nozzle  gross  thrust  coefficient  decreases  significantly.  A  factor  which 
makes  this  transonic  dip  all  the  more  important  is  the  corresponding  increase 

in  the  transonic  drag  coefficient.   Typical  data  for  a  plug  nozzle  are  shown  in 

(18) 
Figure  "*6,  which  is  copied  from  the  NASA  Memorandum  by  Harrington.      Part  (a) 

is  of  interest  here  while  part  (b)  is  not.   Two  flow  phenomena  cause  the  dip 

in  Cp  .   One  is  the  change  in  plug  thrust.   Figure  36(d)  shows  the  behavior 

of  plug  thrust.   Another  influence  is  boattail  drag  identified  in  Figure  36(c) 

as  primary  flap  drag.   Plug  thrust  decreases  and  boattail  drag  increases  in 

the  transonic  region. 

Is  there  anything  that  can  be  done  to  correct  the  loss  of  the  transonic 

V 

(19) 
Rabone     shows  that  by  going  to  small  plug  angles,  some  of  the  loss  of 

plug  thrust  can  be  recovered,  although  there  is  still  a  dip  near  M  =  1.2. 

(18) 
Harrington     uses  a  translating  shroud  that  has  a  pronounced  influence  on 

the  boattail  drag;  however,  the  loss  of  a  plug  thrust  partially  remains. 

An  alternate  approach  is  to  burn  externally  on  the  boattail.   The  direct 
thrust  produced  by  high  pressure  on  the  boattail  has  very  poor  SFC.   However, 
a  flow  interaction  may  occur  which  greatly  increases  thrust.   This  phenomenon 
has  been  termed  "wave  trapping"  by  Fuhs.   Figure  37  illustrates  some  essential 
features  of  the  interaction. 

With  heat  addition  on  the  boattail,  the  turning  angle,  9,  of  the  primary 
nozzle  flow  is  less.   See  the  angle  0  in  the  region  labelled  (l)  in  Figure  37 . 
At  the  slipstream  between  the  external  flow  and  the  nozzle  flow,  the  expansion 
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Nozzle  gross  thrust  coefficient,  C 
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waves  are  reflected  as  expansion  waves  without  heat  addition.   These  reflected 
expansion  waves  further  reduce  the  pressure  on  the  plug.   The  reflection  is 
indicated  at  point  (2)  in  Figure  T?.   With  heat  addition,  the  expansion  wave 
is  reflected  as  a  compression  wave  increasing  pressure  on  the  plug. 

In  region  ("*)   without  heat  addition  the  pressure  is  low,  and  the  Mach 
number  is  large.   The  oblique  shock  {/+)   has  a  shallow  angle.   It  is  necessary 
to  match  slipstream  pressure  and  flow  direction  for  both  the  external  and 
nozzle  flows.   At  the  match  condition,  the  pressure  rise  across  the  obi ique 
shock  is  small,  giving  large  expansion  of  the  nozzle  primary  flow. 

With  heat  addition  there  are  several  changes  to  the  flow  in  region  (3). 
Pressure  is  increased.   The  flow  is  deflected  to  larger  values  of  9.   The  Mach 
number  is  lowered  in  region  (3)  due  to  heat  addition.   Either  a  weak  or  a  strong 
solution  of  the  oblique  shock  occurs.   When  the  strong  shock  solution  occurs, 
the  external  flow  downstream  of  the  oblique  shock  is  high  subsonic.   Expansion 
waves  reflected  off  the  plug  surface  are  not  transmitted  into  the  external 
flow.   Hence  waves  are  trapped  and  are  reflected  as  compression  waves.   The 
compression  waves  maintain  a  high  plug  pressure.   Having  looked  at  the  inter- 
action qualitatively,  let  us  now  examine  the  results  of  a  sample  quantitative 
calculation. 

W-.ve  Trapping;  Quantitative  Example 

(IS) 
The  same  geometry  plug  nozzle  was  selected  as  that  tested  by  Harrington. 

The  primary  flap  forming  the  boattail  was  turned  inward  at  an  angle  of  l7  .   A 

freestream  Mach  number  of  1.31  was  chosen  for  the  calculations;  one  reason  for 

selecting  this  value  of  Mach  number  was  to  keep  the  external  flow  supersonic 

throughout  when  there  was  no  heat  addition.   A  nozzle  pressure  ratio  of  6  was 
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chosen.   According  to  Figure  5,  at  maximum  thrust  a  higher  NPR  would  be  more 

appropriate  for  Kn  =  1.31.   NPR  of  6  is  characteristic  of  cruise.   In  order  that 

standard  gas  tables  could  be  used,  the  ratio  of  heat  capacities,  Y,  was  set 

equal  to  1.4.0. 

To  obtain  the  flow  field,  a  planar,  two-dimensional,  finite,  wave  calcula- 

(15) 
tion  procedure  outlined  in  Section  12.9  of  Liepmann  and  Roshko     was  used. 

When  the  zone  of  influence,  as  defined  by  Oswatitsch     in  Section  3.25  of  his 

book,  spans  a  distance  small  compared  to  flow  radius,  the  two-dimensional  planar 

procedure  can  give  quite  accurate  answers  for  local  values  of  flow.   Near  the 

plug  nozzle  throat,  this  procedure  should  give  small  errors;  however,  near  the 

plug  tip,  sizeable  error  should  be  expected.   In  addition,  the  change  in  wave 

angle  at  the  intersection  of  finite  waves  was  neglected. 

Results  of  the  calculation  are  shown  in  Figures  38  and  ?9.   Figure  ?8  shows 
the  wave  geometry.   The  dashed  line  is  the  slipstream  between  the  external  and 
the  nozzle  flows. 

Consider  the  E  wave*  separating  regions  (l )  and  (2).   It  is  reflected 

as  an  expansion  wave  from  the  plug,  becoming  an  Eo  wave.   The  Ed  wave  reflects 

from  the  slipstream  immediately  downstream  of  region  (4).   The  value  of  the 

reflection  coefficient  is  ♦  0.3.   This  means  the  reflected  wave  is  E  with  a 

r 

strength  of  0.3  of  that  of  the  incident  Ep  wave.   At  the  next  reflection  from  the 
slipstream,  the  coefficient  is  «■  0.12,  which  means  continued  expansion.   At  the 
next  reflection,  the  coefficient  is  -  0.08;  there  is  an  extremely  weak  com- 
pression wave  reflected  from  the  slipstream.   As  a  consequence  of  the  expansion, 


*E  means  expansion,  and  C  means  a  compression  wave.   Subscript  r  is  for  a 
right  running  wave,  whereas  subscript  I   is  for  a  £eft  running  wave. 
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the  pressure  along  the  plug  drops  below  ambient  as  shown  in  Figure  39,   Pressure 
below  ambient  on  the  plug  results  in  a  negative  plug  thrust.   Integration  of  the 
curve  shown  in  Figure  39  yields  a  value  of  the  ratio,  plug  thrust  to  ideal  nozzle 
thrust,  equal  to  -  0.09"7.   Direct  comparison  of  this  calculation  with  Harrington's 
data  is  not  possible  since  all  his  tests  at  M  =  1.31  were  run  with  a  shroud 
extending  partially  over  the  plug.   However,  the  value  of  plug  thrust  has  the 
correct  algebraic  sign  and  the  correct  order  of  magnitude.   In  fact,  it  is 
better  than  the  correct  order  of  magnitude;  it  is  probably  within  50  per  cent 
of  being  correct. 

Now  let's  add  heat  to  the  region  near  the  boattail  as  shown  in  Figures  317 
and  40.   The  length  of  the  boattail  is  L.   The  thickness  of  the  heat  addition 
region  is  I.      The  pressure  increase  due  to  heat  addition  is  obtained  from 
Equation  (41).   Evaluation  of  the  integral  yields 

.  (T  -l)Mht  (   fi) 

r        ap 

Equation  (136)  can  be  solved  for  the  heat  addition  h,  BTU/ volume  second,  and 

then  multiplied  by  I/t. 

Tj>,     aEp'L       BTU ,_««* 

"h=  (Y  -l)M  (length)(sec)  (l?) 

Multiplication  by  ^r,  i.e.,  the  circumference  of  the  boattail,  leads  to  the 

time  rate  at  which  heat  is  released.   Let  H  be  the  heating  value  of  the  fuel 

and  mf,  the  fuel  flow  rate.   The  heat  release  rate  is  given  by 

V "  TT%  (1,8) 

The  specific  fuel  consumption,  SFC,  is  given  by 

3600m 
5FC  -   -g-t  (1*» 
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FIGURE  40.  HEAT  ADDITION  REGION  AT  BOATTAIL. 
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when  AD  is   the   change   in   boattail  drag.      The   change   in  boattail  drag    is 
AD  =    Lp1    sin  o;    the   angle  a    is   defined   in  Figure  4-0.      Combining  Equations 
(138)   and    (i"3?)   along  with  the  expression  for  AD  gives 

SFC  -    ,        *°ffi    .  (140) 

(Y  -  1)MH  sin  o 

Inserting  numerical  values,  a  =  1000  ft/sec,  p  =  "yM  -  1  =  1.615,  Y  =  1.4, 
M  =  1.92,  H  =  20,000  BTU/lb  ,  and  a  =  17°,  and  dividing  by  778  ft-lbf/BTU 
leads  to  SFC  =  1.685  per  hour.   This  is  not  particularly  an  exciting  value  for 
SFC;  however,  for  the  strong  solution  the  heat  addition  causes  wave  trapping  and 
significantly  increases  plug  thrust.   Note  that  SFC  is  not  dependent  on  the 
amount  of  heat  added.   Less  heat  gives  less  thrust;  more  heat  gives  propor- 
tionately more  thrust. 

Now  let's  look  at  the  flow  on  the  plug.   To  start  the  solution,  the  slip- 
stream angle,  9,  must  be  determined.   Figure  41  is  a  plot  of  the  curves  showing 
pressure  in  the  external  flow  as  a  function  of  9  and  the  pressure  of  the  nozzle 
flow  also  as  a  function  of  9.   The  ordinate  is  pressure  behind  the  oblique 
shock  (or  downstream  of  the  expansion  fan  for  the  nozzle  flow)  divided  by 
freestream  stagnation  pressure.   There  are  two  possible  solutions;  the  strong 
solution  was  chosen  for  this  example.   Values  are  9  =  +0.3  ,  subsonic  external 
flow,  and  a  shock  wave  angle  of  77  .   Knowing  9  at  the  boattail-slipstream 
junction  permits  one  to  start  the  solution  using  finite  waves. 

The  wave  geometry  along  the  plug  is  shown  in  Figure  42.   Since  the 
reflection  coefficients  at  the  slipstream  have  the  value  -  1.0,  the  flow 
becomes  periodic.   The  pressure  along  the  plug,  which  is  shown  in  Figure  43,  is 
everywhere  larger  than  ambient  pressure.   Integration  of  the  pressure  gives 
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plug  thrust      =  o  22? 
ideal  nozzle  thrust 

If  heat  is  added  so  that  p'/Prp/->  =  .027  on  the  boattail,  the  value  of  SFC 

becomes  less  than  0.1  per  hour  .'   In  Equation  (l?9),  AD  would  now  include 

both  the  decrease  in  boattail  drag  and  the  increase  in  plug  thrust. 

Wave  Trapping ;  Discussion  of  Numerical  Example 

In  the  course  of  making  a  numerical  experiment,  there  are  numerous 
assumptions  and  decisions  to  make.   Having  completed  this  part  of  an  analysis 
of  the  influence  of  heat  addition,  one  has  new  perspective.   The  decrease  in 
boattail  drag  and  the  corresponding  SFC  are  fairly  straightforward  calcula- 
tions, and  the  value  of  SFC  =  1.685  is  probably  fairly  accurate. 

Considering  the  plug  thrust  increment  due  to  heat  addition,  it  is  apparent 
that  a  crossroads  occurs  at  the  point  of  taking  a  weak  or  strong  solution. 
See,  once  again,  Figure  41.   For  this  example,  the  strong  solution  was  taken. 

The  strong  solution  gives  a  subsonic  external  flow  at  least  initially. 
That  fact,  of  course,  is  the  basis  for  wave  trapping.   Somewhere  along  the  plug, 
the  external  flow  may  accelerate  from  subsonic  flow  to  supersonic  flow.   Once 
the  external  flow  becomes  supersonic,  the  waves  are  no  longer  trapped  in  the 
nozzle  flow.   Determination  of  the  subsonic  flow  and  the  change  to  supersonic 
flow  is  extremely  difficult. 

If  one  had  taken  the  weak  solution  instead  of  the  strong  solution,  the 
pressure  distribution  along  the  plug,  at  least  according  to  current  estimates, 
would  be  qualitatively  similar  to  Figure  ?9. 

Changing  NPR  raises  the  curve  labelled  "NOZZLE  FLOW"  in  Figure  41. 
As  the  nozzle-flow  curve  is  raised,  the  weak  and  strong  solutions  converge  to  a 
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single  solution  at  the  point  of  tangency.   See  Figure  4A(a)  and  4-4.  (b). 
This  single  solution  is  a  strong  shock  solution.   Increasing  NPR  beyond  the 
point  of  tangency  raises  the  curve  for  "NOZZLE  FLOW"  so  that  an  intersection 
of  the  nozzle  flow  curve  and  external  flow  curve  does  not  occur.   The  oblique 
shock  probably  moves  upstream  to  the  point  where  the  boattail  is  formed. 
This  case  is  illustrated  in  Figure  4-4- (c ).   This  would  give  high  pressure  on 
the  boattail  but  probably  would  give  low  plug  thrust. 

Adding  heat  on  the  boattail  has  two  influences:   (l)  Due  to  change  in  9, 
the  curve  for  external  flow  in  Figure  4-1  is  shifted  to  the  right;  and  (2)  due 
to  increased  pressure  ahead  of  the  shock  wave,  the  curve  is  moved  vertically. 
This  is  illustrated  in  Figure  45  (a).   For  large  NPR  and  small  heat  addition, 
there  may  be  a  solution  as  shown  in  Figure  4.5(c).   Additional  heat  release 
causes  a  shift  in  the  curve  making  possible  both  weak  and  strong  solutions. 
See  Figure  4.5(d). 

The  validity  of  the  flow  shown  in  Figures  4-2  and  4.3  is  doubtful  due  to 
the  fact  that  the  external  flow  probably  accelerates  from  subsonic  to  super- 
sonic.  The  strong  solution  may  not  occur.   Consequently,  one  is  skeptical 
of  the  SFC  =  0.1.   Additional  analysis  is  necessary  to  obtain  the  plug 
pressure  distribution  with  heat  addition.   Both  the  weak  solution  and  refine- 
ment of  the  strong  solution  should  be  considered.   One  would  expect  the  heat 
addition  to  have  a  favorable  influence  on  plug  thrust. 
Base  Flow  Problem  with  Heat  Addition;  Flanar  Flow 

One  application  of  external  burning  which  has  been  investigated  both 
experimentally  and  theoretically  is  the  modification  of  base  pressure  by  heat 
addition.   The  heat  may  be  added  in  the  base  recirculation  zone,  in  the  viscous 
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FIGURE  44. 

INFLUENCE  OF  CHANGING  NOZZLE  PRESSURE  RATIO  WITH  FIXED  HEAT  ADDITION. 
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shear  layer,  or  in  the  inviscid  flow  adjacent  to  the  shear  layer.   In  1°67, 

f  1  & ) 

Billig  "    prepared  an  excellent  summary  of  work  to  date.   More  recently, 

(20) 
Roberts     has  written  a  survey  on  the  subject. 

It  is  worthwhile  to  examine  the  base  flow  from  an  integral  point  of  view. 

This  is  illustrated  in  Figure  4.6.   The  control  volume  for  the  momentum  theorem 

is  shown  as  dashed  lines  with  various  surfaces  identified  by  a,  b,  c,  .  .  .  . 

(15) 

Using  notation  of  Chapter  7  of  Liepmann  and  Roshko,     the  momentum  theorem  is 


I  pu.  ( d.n.  )dA  +    pn. 

J    -  JJ         J    x 


dA  =  0 


(HI) 


One  can  add  an  integral  of  ambient  pressure  over  the  control  volume.   For 
base  pressure  i  =  1. 


pu,(u1n  ♦  u-n  )dA  + 


(p  -  p^n-jdA  =  0 


(K2) 


For  the  control  volume  shown  in  Figure  4-6 ,  there  are  four  integrals. 


J 


(pb  -PJ(-l)dA  + 


r 


pu  (-  l)dA  + 


be 


(p  -  P^MA  +     pi^dA  =  0   (K3) 
cd  de 


(pb  -  Poe)dA  =  pi^dA  ♦  -  pu^dA  ♦        \    (p  -  P00)n1dA  (1U) 

ab  de  ^-^.  be       ^-^  cd  /^-^. 

Cl)  ©  Cl) 

Heat  addition  in  the  inviscid  flow  above  the  edge  of  the  shear  layer  will  alter 
the  shape  of  the  streamline  from  c  to  d.   Changing  streamline  cd  will  not 
influence  integral  (  2 J  above.   It  is  not  readily  apparent  how  integral f  l) 
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changes  due  to  modification  of  streamline  cd.   Of  course  integral  (  1j   is 
directly  influenced. 

To  gain  some  insight  into  the  anticipated  change  in  integral  (  lj,    one 
can  use  the  continuity  equation  which  is 

pu-jdA  =   1  pt^dA  (H5) 

be         de 
Since  integral  over  surface  be  doesn't  change  when  streamline  cd  is  changed, 

this  means  the  integral  over  de  does  not  change.   Any  variation  in  integral  (  lj 

2 

of  Equation  (144)  is  due  to  variations  in  p  and  u,  ;  however,  Equation  (14-5) 

tells  us  any  change  in  p  is  balanced  by  a  compensating  change  in  u. . 

If  the  flow  above  the  edge  of  the  shear  layer  illustrated  in  Figure  4-6 
is  inviscid  and  if  only  one  family  of  characteristics  is  significant,  then 
one  can  relate  local  streamline  slope  to  pressure.   To  demonstrate  the  influence 
of  the  shape  of  streamline  cd  on  base  pressure,  a  sample  problem  was  worked. 
This  is  illustrated  by  Figure  4-7  .   Streamline  cd  is  used  with  the  same  lateral 
displacement  but  with  large  turning  angle  in  the  upper  diagram.   The  pressure 
distribution  along  the  streamline  is  shown  in  the  graphs  below  each  streamline. 
From  the  curves  one  can  evaluate  the  integral  or  equivalent  series 


I  — n  dA  =  } 

cd  j 


-*W.  (H6) 


where  AA.  is  the  element  of  area  projected  in  the  x,  direction.   Numerical 

J 

values  for  I  are  12.8  for  upper  curve  and  21.7  for  lower  curve.   Relating 
Equation  (14.6)  to  Equation  (144)  gives 
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FIGURE  47    PRESSURE  ON  EDGE  OF  SHEAR  LAYER 
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oo  . 


(-  -  -^)dA  = 


PT   PT 


(p  -  pJryiA  =  pT(l  -  52.8)  (U7) 

cd  cd 

The  sharp  turning  curve  gives  a  value  of  -  40pT  for  integral  (  3  J  of 
Equation  (144),  whereas  the  gradual  turning  curve  gives  -  31. lp^.   The  base 
pressure  would  be  greater  with  the  gradual  turning  streamline  cd.   It  is 
obvious  from  this  analysis  that  one  wants  to  maintain  a  high  pressure  on 
streamline  cd  to  increase  base  pressure.   The  heat  addition  zone  shown  in 
Figure  48  will  do  that. 

Now  let's  return  to  the  problem  of  how  much  pressure  increase  results 
from  the  heat  addition  zone.   Consider  a  heat  addition  zone  located  above  a 
slipstream  as  shown  in  Figure  49.   This  is  a  model  to  represent  the  flow 
depicted  in  Figure  48.   One  replaces  the  lateral  gradient  in  Mach  number 
with  a  slipstream  having  supersonic  flow  above  and  subsonic  flow  below.   With 
a  lateral  Mach  number  gradient,  as  in  Figure  48,  the  characteristics  reflect 
from  the  sonic  line  in  the  form  of  a  cusp. 

A  paragraph  will  now  be  used  to  explain  the  notation  E  ,  C  ,  Eb,  Cb. 
The  symbols  E  and  C  are  for  Expansion  and  Compression  waves.   The  subscripts 
r  and  I   represent  right  running  and  left  running  waves  respectively. 

In  Figure  49  the  heat  addition  zone  has  a  leading  edge  swept  at  angle  ji. 
This  gives  a  finite  wave  C  .   The  pressure  along  the  slipstream  remains  at 
pressure  p  .   The  waves  from  the  heat  addition  zone  and  the  reflected  waves 
turn  the  flow  while  keeping  pressure  constant.   One  way  to  describe  this  is  to 
say  that  the  compression  waves  from  the  heat  addition  zone  have  cancelled  the 
expansion  waves  which  would  have  originated  due  to  kink  in  streamline  at 
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FIGURE    48. 
INTERACTION  OF  HEAT  ADDITION  ZONES  WITH   BASE  FLOW. 
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int  (A).   The  turning  angle  is  2©  at  point  ( A)   or  twice  the  turning  angle  due 


po: 

to  C  .   The  reflection  has  effectively  doubled  the  influence  of  the  waves  from 

r 

the  heat  addition  zone. 

The  doubling  of  effectiveness  can  be  seen  in  the  following  analysis:   In 
free  unbounded  space  the  turning  angle  is  given  by 

tan  q  _  e  =   v'  =  zi  =   (r  -  i)m>  ,  g) 

tan  y  -  y  -  n  +  u,   0     ^pU  U4«; 

Equation  (14-8)  is  obtained  from  geometry  and  Equation  (4-0).   From  Equation  (4-1) 

p,  .  (Y  -al)Mhb  {U9) 

where  symbol  b  is  defined  by  Figure  49.   Due  to  turning  of  flow  by  an 
amount  20,  the  increment  of  pressure  decrease  would  be 

p*  =  .  2n^9  (150) 

From  the  flow  field  analysis  illustrated  in  Figure  4-9,  the  following 
equation  must  be  true 

2p«  +  prt  =  0  (151) 

Combining  Equations  (14-8),  (14-9),  (150 ),  and  (151;,  one  can  verify  that 
Equation  (151 )  is  in  fact  true. 

The  mechanism  for  altering  base  pressure  by  heat  addition  in  the  adjacent 
inviscid  flow  is  apparent  if  one  retraces  the  discussion  associated  with 
Figures  4-6,  4-7,  and  4-8.   To  summarize,  the  streamline  cd  must  turn  inward 
toward  the  axis  or  plane  of  symmetry.   Without  heat  addition,  this  turning 
causes  the  pressure  to  drop.   Sharp  turning  gives  low  pressure  and  low  base 

pressure;  see  Equation  (144-)  and  Figure  4-7«   Heat  addition  allows  this  inward 

(21) 
turning  without  loss  of  pressure  as  shown  in  Figures  48  and  4-9.   Strahle's 

analysis  suggests  a  new  way  to  overcome  base  drag. 
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Base  Flow  Froblem  with  Heat  Addition;  Axisymmetric 

(p2)  (2^) 

The  analysis  of  Fein     and  Strahle     considered  a  two-dimensional 

planar  geometry,  i.e.,  wing-like  objects.   Frojectiles  are,  of  course, 

axisymmetric.   This  section  examines  some  of  the  changes  to  be  expected  when 

the  planar  results  are  replaced  by  calculations  for  the  axisymmetric  case. 

(15) 
The  method  outlined  in  Liepmann  and  Roshko     was  used  to  calculate  the 

flow  field  for  axisymmetric  heat  addition.   To  verify  the  procedure  for 

handling  the  characteristics  near  the  axis,  a  model  of  a  spherical,  radially 

expanding  flow  was  used.   The  results  of  that  calculation  are  shown  in  Figure  50. 

Since  the  flow  is  known,  one  knows  Mach  number  at  any  radius  R.   The  Mach 

number  calculated  from  the  characteristics  solution  is  shown  in  the  third 

column  cf  the  data  summary  in  Figure  50.   The  Mach  number  based  on  the  value 

of  R  and  flow  area  is  shown  in  the  first  column. 

A  planar  radial  flow  has  an  area  variation  linear  in  R,  whereas  the 

2 

axisymmetric  radial  flow  has  a  R  dependence  for  area.   This  can  be  represented 

by 

M  _   AR  AA  _    ^AR  ,nco> 

T~T  T-  *R  (152) 

Planar  Axisymmetric 

Partial   solutions  of  the  flow  field  with  heat  addition  were  obtained  for 

two  cases.      Both  cases  had  an  annular  heat  addition  zone  as   illustrated   in  the 

lower  left-hand   corner  of  Figure  51.      The  heat  addition  was   sufficient  to   cause 

a  maximum  turning  angle  of  0         =    -  6     at  a  radius   of  6   inches.      In  the  first  case, 
&        6  max 

the   heat  zone  was   bounded  by  6  ^  r  ^  ^  an^  M  z  O-      ^n  the  second   case,    it 
was  6  ^  r  ^  8  and  0  ^  z   ^3.      For  the  first  case,    shown  in  Figure  51,   Mq  =    -\]2. 
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FIGURE   50.    CHARACTERISTICS  SOLUTION   FOR  AXISYMMETRIC 
RADIAL  EXPANSION. 
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For  the  second  case,  M~  was  chosen  as  2,  and  the  radial  extent  of  the  heat 
addition  zone  was  increased  to  spread  out  the  waves.   This  was  done  to  avoid 
merging  of  the  characteristics  to  form  a  shock  wave. 

Pressure  versus  radius  is  plotted  in  Figure  51-   The  mid-characteristic, 
which  has  points  1,  2,  and  3,  has  an  imperceptible  increase  in  pressure  as 

radius  decreases.   The  characteristics  converge  to  form  a  shnck  wave.   When  this 

( 23  )  (2A) 

happens,  the  techniques  outlined  in  the  book  '    and  NACA  TN     by  Ferri  are 

applied. 

Along  the  characteristic,  which  starts  at  r  =  6"  with  9    =  -  6  ,  the 

max       ' 

pressure  rises  slightly  to  a  radius  of  r  =  2".   Between  0  ^  r  ^  2,  the  pressure 
increases  rapidly.   Due  to  the  large  mesh  size,  the  pressure  variation  at 
points  12,  11,  and  14  may  not  be  correct.   Certainly  the  end  point,  number  14., 
is  correct.   The  dashed  line  represents  the  curve  that  would  probably 

obtained  with  finer  mesh  size  and  a  more  accurate  treatment  of  the  shock  waves. 

2 

Note  that  for  M„  =  2,  even  a  normal  shock  wave  is  nearly  isentropic,  which 

permits  treatment  of  shock  waves  as  finite  isentropic  compression  waves. 
The  conclusions  from  Figure  51  are  as  follows: 

(a)  In  the  axisymmetric  geometry,  compression  waves  tend  to  merge  more 
quickly. 

(b)  The  increase  of  pressure  with  decreasing  radius  is  very  small  until 
the  axis  is  approached. 

Now  let's  look  at  the  second  example  which  was  chosen  to  avoid  a  shock  wave. 

The  characteristics  for  the  second  case  are  shown  in  Figure  52.   Increasing 
the  spacing  of  the  waves  did  not  completely  avoid  a  shock  wave.   A  normal  shock 
extends  from  points  1?  and  1U   to  the  axis. 
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The  pressure  as  a  function  of  radius  is  shown  in  Figure  53  for  the  mid-wave 
(wave  starting  at  r  =  6  with  6  =  -  3  )  and  the  terminal  wave  (wave  starting  at 
r  -  6  with  8  =  -  6  ).   As  before  in  the  previous  example,  the  pressure  does  not 
increase  significantly  until  very  close  to  the  axis.   This  is  true  for  both  the 
mid-wave  and  the  terminal  wave.   The  pressure  is  shown  as  a  dash-dot-dash  line 
between  point  1/+  and  the  axis  since  this  was  not.  calculated  in  detail. 

Figure  5U   gives  information  about  the  streamlines  and  flow  deflection. 
Right  on  the  axis,  of  course,  the  streamline  is  the  axis.   There  is  large 
turning  near  the  axis  giving  small  radius  of  curvature  for  the  streamline. 
This  balances  the  large  radial  pressure  gradient  near  the  axis. 

The  conclusions  from  the  second  case  are  the  same  as  the  first  case. 
Preliminary  Comments  Based  on  a  Cursory  Look  at  Axisymmetric  Case 

The  change  in  base  drag  is  given  by 

ADfe  =  ApbAb  (151) 

Each  of  the  two  terms,  Ap,  and  A,  ,  will  now  be  discussed.   The  results  shown 

b      b 

in  Figures  $1,  52,  and  53  suggest  that  Ap,  may  not  be  increased  significantly 

as  a  direct  result  of  changing  geometry  from  planar  to  axisymmetric.   However, 

the  results  shown  in  Figures  51,  52,  and  54-  more  closely  apply  to  pressure  at 

the  dividing  streamline  and  not  to  base  pressure.   The  flow  in  the  base  region 

may  magnify  the  pressure,  just  as  the  initial  pressure  due  to  heat  addition  at 

r  =  6  is  increased  slightly  due  to  axisymmetric  geometry. 

Compression  originating  with  the  heat  addition  will  be  reflected  at  the 

sonic  line.   See  Figure  4-8  and  point  (\k)   of  Figure  55.   These  compression  waves 

will  be  reflected  from  the  wake  of  the  heat  addition  zone;  see  point  Qi)   of 

Figure  55.   The  reflection  may  be  negative,  i.e.,  changing  the  compression  waves 

to  expansion  waves.   This  negative  reflection  would  not  be  favorable. 
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FIGURE  53.  PRESSURE  DISTRIBUTION  ALONG  CHARACTERISTICS  OF  FIGURE  52. 
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(B)  FLOW  DEFLECTION  AS  A  FUNCTION  OF  RADIUS 


FIGURE  54.  STREAMLINES  AND  FLOW  DEFLECTION. 
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Now  let's  look  at  the  area  term.   Consider  the  same  amount  of  heat,  Q, 

released  in  both  planar  and  axisymmetric  case.   In  the  axisymmetric  case,  Q  is 

distributed  along  a  circle  of  radius  r~  giving  A/2^r„,  BTU/length  sec.   In  the 

planar  case,  Q  is  distributed  along  two  lines  of  length  "r,.   One  line  is  at  the 

top,  and  the  other  at  the  bottom.   See  Figure  56.   For  the  planar  case,  the 

heat  release  is  Q/^r^,  BTU/length  sec,  which  is  identical  to  the  axisymmetric 

case.   The  area  influenced  by  the  heat  release  is  in  the  ratio 

2 
base  area;  axisymmetric  case    "r,     r 

base  area;  planar  case      ^r.r    2r~ 

Making  r, /r~  small  helps  increase  pressure  at  r,  due  to  heat  addition;  see 
Figure  56.   Making  r, /r  large  makes  the  base  area  very  small  in  axisymmetric 
case  as  compared  to  the  planar  case.   One  factor  helps;  the  other  hurts. 

As  a  final  remark,  it  should  be  stated  again  that  these  are  preliminary 
thoughts  and  conclusions. 
Spin  Recovery  Wing  Heat  Addition 

Some  modern  fighter  and  attack  aircraft  can  be  flown  into  spins  of 
such  nature  that  special  equipment  is  needed  for  recovery.   The  drag 
parachute  may  be  deployed  to  lift  the  aircraft  tail.   Before  normal  flight 
is  restored,  the  parachute  must  be  jettisoned.   An  alternate  method  of 
obtaining  a  torque  on  the  aircraft  would  be  heat  addition.   Figure  5?  illustrates 
a  modern  fighter  in  a  spin.   The  engines  run  normally  and  produce  thrust.   By 
spraying  fuel  into  the  exhaust  stream,  the  heat  addition  deflects  the  jet  down- 
ward.  A  torque  is  developed  which  tends  to  lift  the  tail.   Let's  compare  a 
parachute  recovery  with  heat  addition. 
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FIGURE  56.  COMPARISON  OF  BASE  AREA  FOR  PLANAR  AND  AXISYMMETRIC  GEOMETRY. 
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FIGURE    57.   SPIN  RECOVERY  USING  HEAT  ADDITION 
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6     2 
Assume  the  aircraft  has  a  moment  of  inertia  of  6  x  10  lb  ft  about  the 

pitch  axis.   Assume  the  installed  thrust  is  35,000  lb  with  afterburner.   If 

the  flow  can  be  deflected  by  6  with  heat  addition,  the  force  tending  to  lift 

the  tail  is  4.000  lb.   If  the  aircraft  has  a  moment  arm  of  22  feet,  the 

restoring  torque  is  8.8  x  10  lb  ft.   Time  required  to  point  the  aircraft 

downward,  i.e.,  turn  90  in  pitch,  is 

t  =  -y  Y"  =  ^  sec 
If  the  sink  rate  is  1^0  ft/sec,  the  recovery  altitude  is  64.00  ft. 

One  can  estimate  the  amount  of  heat  to  be  added  for  a  deflection  of  6 
using  Equations  (119)  through  (129).   The  temperature  of  exhaust  must  be 
increased  by  30  per  cent. 

Now  consider  a  parachute  of  10  ft  diameter  with  a  drag  coefficient  of  2. 
The  drag  at  a  sink  rate  of  140  ft/sec  will  be 

D  =  ^fp^AC.  =  ^200  lb 

This  is  comparable  to  the  force  due  to  thrust  deflection.   Performance  of 
parachute  and  heat  addition  is  comparable. 

SUMMARY  AND  CONCLUDING  DISCUSSION 
Use  of  an  area  rule  or  quasi  area  rule  implies  optimization.   The  area 
rule  can  be  used  to  translate  an  optimum  axisymmetric  body,  e.g.,  von  Karman 
ogive,  to  a  three-dimensional  body.   The  optimum  axisymmetric  body  shape  is  found 
by  application  of  the  calculus  of  variations  to  the  source  distribution  function, 
f^x).   When  there  is  heat  addition,  there  is  a  new  function  to  optimize,  Q(x,r); 
see  Figure  55.   Both  f (x)  and  Q(x,r)  must  be  simultaneously  optimized.   The 
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ax i symmetric,  small  perturbation,  solution  for  heat  addition  has  never  been 
obtained;  it  is  an  essential  building  block  of  the  quasi  area  rule. 

When  the  body  develops  lift  with  aerodynamic  circulation,  there  is  another 
complexity  to  the  analysis.   Vorticity  due  to  lift  must  be  included  in  the  analysis. 
There  may  be  an  optimum  balance  between  lift  generated  aerodynamically  and  lift 
generated  by  heat  addition. 

To  apply  the  optimization  technique  to  turbojets,  ramjets,  or  other 
internally  burning,  air  breathing,  engines,  a  different  analytical  technique 
is  required.   A  combined  thrust-drag  technique  was  briefly  introduced.   This 
method  uses  an  energy  disc. 

By  discussing  both  one-dimensional  and  two-dimensional,  planar,  heat 
addition,  the  connection  between  heat  fronts,  combustion  fans,  and  one-dimensional 
heat  addition  was  demonstrated.   Heat  addition  in  the  planar  case  may  be  represented 
either  by  line  sources  or  by  a  volume  distribution  function,  Q(x,r).   In  several 
ways  heat  addition  is  similar  to  the  flow  generated  by  a  solid  body.   The 
pressure  en  the  flat  plate  of  Figure  32  is  the  same  as  that  caused  by  a  wedge 
of  angle  9.   Pressure  and  velocity  perturbations  are  identical  for  both  solid 
body  and  heat  addition.   Differences  exist  between  flow  over  solid  bodies 
and  flow  through  a  region  of  heat  addition.   One  difference  is  drag. 

There  are  three  approaches  to  calculate  thrust  minus  drag  in  an  inviscid 
flow:   integration  of  pressure  over  the  body  surface,  momentum  control  volume, 
and  wave  drag.   Integration  of  pressure  is  the  most  direct  of  the  methods. 
Figures  22  and  27  illustrate  the  wave  energy  and  momentum  approaches  for  a 
thrusting,  lifting,  planar  airfoil.   Figures  20,  23,  and  26  illustrate  the  wave 
and  momentum  approaches  for  an  isolated  heat  addition  zone  in  an  infinite  medium. 
The  merit  of  a  particular  method  depends  on  the  geometry. 

128 


Applications  of  external  burning  include  lifting,  thrusting  airfoils,  base 
or  boattail  pressure  modification  for  an  exhaust,  development  of  base  thrust 
for  projectiles,  and  spin  recovery.   External  burning  for  airfoils  and  projectiles 
has  been  studied  extensively  both  theoretically  and  experimentally.   A  new 
phenomenon  occurs  for  the  exhaust  problem;  this  is  wave  trapping.   A 
preliminary  analysis  indicates  heat  addition  for  spin  recovery  is  comparable 
to  use  of  a  parachute. 
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